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Abstract 



The theory of generalized local scale invariance of strongly anisotropic scale invariant systems 

S proposed some time ago by Henkel [Nucl. Phys. B 641, 405 (2002)] is examined. The case of 

so-called type-I systems is considered. This was conjectured to be realized by systems at m-axial 
Lifshitz points; in support of this claim, scaling functions of two-point cumulants at the uniaxial 
Lifshitz point of the three-dimensional ANNNI model were predicted on the basis of this theory 
and found to be in excellent agreement with Monte Carlo results [Phys. Rev. Lett. 87, 125702 
Pi \ (2001)]. The consequences of the conjectured invariance equations are investigated, with empha- 

sis put on this case. It is shown that fewer solutions than anticipated by Henkel generally exist 
and contribute to the scaling functions if these equations are assumed to hold for all (positive 
and negative) values of the d-dimensional space (or space time) coordinates (t,r) € M x M d_1 . 
Specifically, a single rather than two independent physically acceptable solutions exists in the 
case relevant for the mentioned fit of Monte Carlo data for the ANNNI model. Renormalization- 



CN 
> 

m . 

r-f\ , group improved perturbation theory in 4 + m/2 — e dimensions is used to determine the scaling 

functions of the order-parameter and energy- density two-point cumulants in momentum space 
to two-loop order. The results are mathematically incompatible with Henkel's predictions ex- 
cept in free-field-theory cases. However, the scaling function of the energy-density cumulant we 
obtain for m — 1 upon extrapolation of our two-loop RG results to d — 3 differs numerically 
little from that of an effective free field theory. 
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1. Introduction 

It is a fact, well established by plenty of experiments and theoretical works, that 
the large-scale physics of systems at critical points can be described by scale-invariant 
continuum field theories [1-3]. In many cases the associated probability distributions 
are besides scale invariant also translation and rotation invariant when expressed in 
appropriate variables. During the past 25 years it has become widely appreciated that in 
those cases where long-range interactions are absent or may be ignored, such invariancc 
under translations, rotations, and scale transformations usually entails the invariancc 
under a larger symmetry group — that of conformal transformations [4-8] . 

The application of conformal invariancc to critical phenomena started long ago with 
a short note by Polyakov [9]. More detailed investigations [10,11] followed soon, but 
the issue did not receive much attention until the seminal work of Bclavin ct al [12,13] 
on two-dimensional conformal field theories. This revealed the enormous potential of 
conformal invariancc, triggering an outburst of research activities, which in turn led to 
impressive success in many other fields, such as bulk, finite-size [14-16], and boundary 
critical behavior [5,6,17-20], polymer physics [21], quantum impurity problems [22], and 
string theories [7,8,23]. 

Conformal transformations locally correspond to combinations of translations, rota- 
tions, and scale transformations with a position-dependent scale factor £(x) that involve 
no shear. Thus, conformal invariance may be viewed as the generalization of a global 
symmetry to a local one. 

There exists a wealth of systems in nature that exhibit global scale invariance of a dis- 
tinct and more general kind, called anisotropic scale invariance (ASI). Its characteristic 
feature is that an anisotropic rescaling of the space separations (or spacetime separations 
in time-dependent phenomena) along various axes by at least two (or more) distinct pow- 
ers of a scale factor £ is required to make such systems statistically self-similar. Familiar 
examples are uniaxial dipolar ferro- and antifcrromagnets at their critical points (see, 
e.g., Refs. [24] and [3, chapter 27.5]), systems at Lifshitz points [25-27], and dynamical 
critical phenomena near and away from thermal equilibrium [28,29], among them driven 



diffusive systems [29] , stochastic surface-growth processes [30] , directed percolation, and 
spreading processes [31,32]. 

In the case of static critical behavior at an m-axial Lifshitz point (LP) in d space 
dimensions, the position vector x — (xi, ■ ■ ■ ,Xd) of Cartesian coordinates x 7 can be 
decomposed as x = (z,r) into the m- and (d — 77i)-dimensional components z = (x a ) 
and r = (xp) with a — 1, . . . , m and j3 — m + 1, . . . , d, respectively. ASI is encoded in 
the transformation property 

O j (l z,lr)=r A *O j (z,r) (1.1) 

of local scaling operators Oj(x) with scaling dimensions Aj, where 9, the anisotropy 
exponent, differs from 1. The obvious analog for time-dependent phenomena involving 
an isotropic rescaling of distances but distinct rescaling of time reads 

Oj(£H,£r) = r A * 3 (t,r) , (1.2) 

where 3 is the so-called dynamic critical exponent. As a consequence of these properties, 
the multi-point correlation functions of such operators take scaling forms. Consider, for 
example, the case of Eq. (1.2), and assume that the systems are translation invariant 
in both time and space, as well as rotation invariant. Together with the presumed ASI, 
these properties imply that the two-point cumulant function of two such operators Oj 
and Ok can be written as 

(Oj(t, r) O k {t>, r')) cum = |r - r'\- A ^ Ak ^ ( j^jl) ' (L3) 

Given the enormous success the use of conformal invariance has had in the study of 
isotropic critical behavior, a natural question to ask is whether global ASI, in conjunction 
with appropriate other global symmetries, such as translation and rotation invariance, 
would again entail more powerful local symmetries that impose useful constraints on the 
scaling functions or even determine them completely. 

This idea has been pursued for many years, in particular, by Henkel who proposed a 
phenomcnological approach termed "local scale invariance (LSI)" in a series of papers 
[16,33,34] and applied it to a variety of systems exhibiting ASI. Postulating a set of 
"axioms of local scale invariance", he suggested that the two-point scaling functions 
fi(u) of various systems exhibiting ASI should satisfy differential equations. According 
to him there should be two classes of local generalizations of ASI: The first, denoted 
type I, should apply to anisotropic scale-invariant equilibrium systems; the second, type 
II, to time-dependent scale-invariant phenomena. As a nontrivial example of type II, 
the relaxational behavior of systems representing the dynamic universality class of the 
so-called stochastic model A [28], following a quench from an initial disordered state to 
the critical point, was suggested. Subsequent analytical calculations based on the e = 
4 — d expansion for model A [35,36] yielded definite, albeit small, violations of Henkel's 
predictions at two-loop order. 

As nontrivial realizations of his type I of generalized ASI, Henkel suggested equilib- 
rium systems at m-axial LP. To check the predictions of his theory, Pleimling and him 
performed extensive Monte Carlo simulations [37] for the two-point correlation function 
of the d = 3 dimensional axial next-nearest-neighbor Ising (ANNNI) model [25-27] at its 
uniaxial LP. In their original application of the phenomenological LSI approach to this 
problem, they assumed that the anisotropy exponent 9 takes its classical value 9 = 1/2. 



Their theory then predicted the scaling function to be a linear combination of two lin- 
early independent solutions of a differential equation, involving a single free parameter 
which they determined from Monte Carlo data for moments. The so-obtained scaling 
function appeared to be in perfect agreement with their Monte Carlo data. However, the 
e expansion about the upper critical dimension d* (to) = 4 + to/2 yields deviations of 9 
from its classical value 1/2 at order e 2 [38-41]. 2 Pade estimates based on these series 
expansions to 0(e 2 ) gave 9 ~ 0.47 for the uniaxial case to = 1 at d = 3. To account for 
such nonclassical values of 9, Pleimling and Hcnkcl [42] generalized their LSI predictions 
for the scaling function by expanding in 9—1/2. They found that the resulting predictions 
remained in agreement with their Monte Carlo data provided a value for 9 sufficiently 
close to 1/2 (0.47 < 9 < 0.5) was chosen. 

Unlike the case of type II, Henkel's predictions obtained via his phenomcnological 
LSI approach have not yet been checked in a systematic fashion by mathematically well 
controlled analytical calculations. The only exceptions we are aware of are mean spherical 
models. Their propagators at the LP are those of massless free field theories. LSI does not 
lead to new nontrivial consequences for them. Hence they are unsuitable for critical checks 
of the predictive power and viability of this approach. In view of the apparent excellent 
agreement of the Monte Carlo data of Ref. [37] with the scaling function obtained by 
the LSI approach we feel that nontrivial checks of this approach through analytical 
calculations for nontrivial models of type-I systems are urgently needed. 

The aim of this paper is to perform such checks. To this end we shall investigate 
standard n-component <fi models for the description of critical behavior at w-axial LP, 
use the e = 4 + to/2 — d expansion to compute appropriate two-point scaling functions, 
and compare the results with the predictions of Henkel's LSI approach. For the sake of 
simplicity, we shall focus our attention in most of our work on the uniaxial case m = 1. We 
begin in Section 2 with a brief review of the main predictions of this theory for the scaling 
functions of such type-I systems. In Section 3 we first discuss the region of validity of the 
suggested invariance equations of the two-point correlation function in position space. 
We then transform these equations to momentum space, and discuss the consequences 
for the scaling forms of the Fourier transformed two-point functions. In Section 4 we 
introduce the standard continuum model representing the universality class of critical 
behavior at to- axial LP in d dimensions. We then investigate the e expansion for the 
energy-density correlation function, using dimensional regularization in conjunction with 
minimal subtraction of ultraviolet (UV) poles. In Section 5.2 we focus on the analytically 
tractable cases d = to + 3 and to = 2. The former has the simplifying feature that 
the scaling function of the free propagator in position space reduces at the LP to a 
Gaussian. This allows us to obtain explicit expressions for the energy-density and order- 
parameter correlation functions to order e 2 and cast them in scaling form. In Section 5.3 
we determine the two-point correlation function of the energy density for the uniaxial 
case to = 1 to two-loop order. A detailed comparison of our results for this and the 
previously mentioned correlation functions at the LP with the predictions of Henkel's 
phenomcnological theory follows in Section 6. It shows that these predictions do not hold 
except in the trivial case of a Gaussian LP. Whenever loop corrections to the correlation 
functions cannot be neglected in the interacting case, the predicted scaling functions 



2 The recently developed large-n expansion for the study of critical behavior at LP, where n is the 
number of components of the order parameter, also gives nonclassical values of 8 for d < d* (m) [40,41] . 



are inconsistent with our findings. Finally, there are four appendices to which we have 
relegated various computational details. 

2. Conjectured properties of scaling functions 

Our objective is to check the predictions of the LSI theory proposed in Rcfs. [16,33,34] 
for strongly anisotropic critical systems of type I. We begin by recalling the basic postu- 
lates on which this theory is based and its conjectured properties of scaling functions. 

To this end, we consider the pair correlation functions of two quasiprimary scaling 
operators Oj(zj, rj), j = 1,2, with zero averages (Oj) and the behavior (1.1) under global 
scale transformations. For the sake of notational simplicity, we focus on the uniaxial case 
m = 1. In order to facilitate comparisons with Hcnkcl's work, we shall follow him and 
denote the one-dimensional equivalent of the variables Zj by tj . We assume translation 
invariance in t space as well as translation and rotation invariance in r space, define the 
scaling dimension 

A=(A 1+ A 2 )/2, (2.1) 

and introduce the notations t = t\ — t 2 , f = i"i — f2, and r = \r\. By analogy with 
Eq. (1.3), the pair correlation functions can be written as 

G jk (t,r) = (0 J (t 1 ,r 1 )O k (t 2 ,r 2 ))=r- 2A n jk (tr- e ). (2.2) 

Their scaling form reflects the invariance under global anisotropic scale transformation 

generated by 

1 2A 

X = -td t --r-d r -—. (2.3) 

a u 

Hence we have 

X Q G jk (t,r) = 0. (2.4) 

To proceed it will be helpful to recall some essentials of Henkel's approach [33,34] 
without going into details. His starting point is the well-known algebra associated with 
Schrodinger invariance. This he generalizes by allowing for values of 7^ 2 and anomalous 
dimensions of scalar quasiprimary fields. He then imposed the requirement that the 
generators yield a finite number of independent conditions when applied to the two-point 
functions of quasiprimary fields. Exploiting the consequences, he was able to identify 
two distinct classes of systems, called type I and type II, respectively. For the type-I 
systems with which we are concerned here, the anisotropy exponent 9 is constrained to 
the fractional values 

9 = N = 2/N, N e N, (2.5) 

so that Eq. (2.4) simplifies to 

X G jk (£, r) = (-td t -^r-d r - NA\ G jk (£, r) = 0. (2.6) 

The other assumptions of Hcnkel are that the Gj k are also annihilated by generators 
denoted as Y"i_]v/2 = (Yi, ■ ■ ■ , Yd-i) and X\ defined via 

Yi-N/2G ik {t,r) =(-td r - ^-rd t N -AG jk (t,r) = (2.7) 

and 



X x G ik {t,r) = (-t 2 d t -Ntr-dr- 2iVAit - a^d^G^r) 

+ 2t 2 X G jk {t, r)+Nr 2 - Y^ N/2 G 3k {t, r) = (2.8) 

where ol\ is a nonzero parameter. 

Note that when N is taken to be an arbitrary real number so that the condition (2.5) 
is not satisfied, the generators Y"i_jv72 an d -Xi involve fractional derivatives. Since defini- 
tions of fractional derivatives d\ other than via Fourier transformation d L t -H- (ik) L are in 
use, the precise definition of these fractional derivatives becomes an issue. Background 
on this matter and Henkel's choice of their definition can be found in reference [34, Ap- 
pendix A] . We shall exclusively have to deal with the above equations in those cases where 
condition (2.5) is satisfied. All derivatives then reduce to conventional partial derivatives 
of first and higher orders. Clearly, any acceptable definition of fractional derivatives d L t 
must reduce to such standard derivatives for nonncgativc integer values of t, i.e., when 
N becomes a natural number. Henkel's choice indeed fulfills this condition. We therefore 
do not have to worry about potential differences resulting from distinct definitions of 
fractional derivatives here and in the following. 

The meaning of the first condition, Eq. (2.6), has already been explained. The second 
condition, Eq. (2.7), reduces in the special cases N = 1 and N — 2 to familiar ones 
implied by invariance under global projective Galilei transformations and rotations, re- 
spectively. The third one, equation (2.8), is reminiscent of the one that follows for systems 
exhibiting isotropic scale invariance in x = {t, r) space from the invariance under Mobius 
transformations. As discussed in Ref. [34, p. 430], the three conditions (2.6)-(2.8) can be 
combined to obtain the constraint 

Ai = A 2 , (2.9) 

unless Gjk = 0. One can therefore put Ai = A and drop the subscript 1 on both Ai and 
a\. 3 Both this constraint and Eq. (2.6) are satisfied by the scaling ansatz 

G jk (t,r)=8 AjAk G(t,r), G(t,r) = r 2A !l«(fr- 2 /*). (2.10) 

Its substitution into Eq. (2.7) then yields a differential equation for the scaling function, 

namely 

rjN-l j \ 

— ^ - v 2 — - (v n {N) (v) = with C - JVA. (2.11) 

dv" l dv J 

Henkcl considers this equation on the interval [0, oo) subject to the boundary conditions 

\im^ N \v) = n (N \0) = J N \ (2.12) 

>( V)(.) « cW^, (2 . 13) 

v— ^oo 

where 0Jq ^0 and Woo are constants. Assuming that N > 2, he arrives at the general 
solutions 

nW( v ) = J2 b i N) v p M v ) ( 2 - 14 ) 

p=0 



3 Reference [34] also uses a parameter 02- However, this is related to ai via «2 = (— 1) N ai in the 
case of type-I systems we are concerned with here. 



with 

where 2-Pat-i is the generalized hypergeometric function, while b p are free parameters. 
Using known theorems [43] about the asymptotic behavior of the functions 2Fn-i{x) 
in the limit x — > 00, he finds that the right-hand side of the solutions (2.14), for general 
values of b p , would diverge asymptotically as 

nW(«) ~ ( va -V N f +1 - NmN - 2 K X p\^(va-V N ) N/{N - 2) ] (2.16) 

v— >oo I 1\ - 

in the large-w limit, where the proportionality constant is a linear combination of the 
coefficients b p . Since such behavior is inconsistent with the boundary condition (2.13), 
he requires that this constant vanishes. This implies the condition 



JV-2 

h( N ) 

r((p+l)/N)T((p + ()/N)\N2J 



V 6™ T{P+1) (^Y /N = (2 17) 



p=0 



l(N) 



which can be used to eliminate the coefficient b N _ 2 . As a consequence, the solutions (2.14) 
become 

N-3 

rt N \v)=^b p N ^Q p N \v), (2.18) 

P =a 

with 

flW M - v v T (v) p\T[(N-l)/N]T[l + (C-2)/N] 

p [) ~ p{] (N-2)\r[(p + i)/N]r[( P + o/N] 

x {a/N 2 ) (p+2 - N)/N v N - 2 F N - 2 {v). (2.19) 

Condition (2.17) ensures the cancellation of the leading exponentially diverging terms 
in the limit v — > +00 of QS N >{v). In order to comply with the boundary condition (2.13), 
no other diverging or non-decaying terms would have to remain in Vv- N >{v) in the limit 
v — >• 00. Provided this is the case, the general solution of equation (2.11) subject to the 
boundary conditions (2.12) and (2.13) involves N — 2 free parameters b p with p = 
0,...,N-3 and is given by Eqs. (2.18) and (2.19). That the boundary condition (2.13) 
is satisfied was confirmed in Ref. [34] by numerical means for N = 4, 5, 6. 

In Appendix B we reconsider in detail the problem of solving Eq. (2.11) subject to the 
boundary conditions (2.12) and (2.13). We prove there the following statements about 
the solutions given by Eqs. (2.18) and (2.19) with general values of the coefficients b p ' : 
When TV = 3,4,5, they comply indeed with the boundary condition (2.13). For general 
integer values N > 7, this boundary condition gets violated by the presence of exponen- 
tially diverging terms in the large-w limit. Requiring the absence of these (subleading) 
divergences imposes further restrictions on the coefficients b p , which reduce their num- 
ber. For example, for N = 7, the general solution of Eq. (2.11) satisfying the boundary 
conditions (2.12) and (2.13) involves only 3 rather than 5 free parameters. The case 
N = 6 is special. As we show in Appendix B, the scaling function fl^ suggested by 



Henkel (for general values of the coefficients b p , p = 0, . . . , 3) violates again the bound- 
ary condition (2.13) but diverges only algebraically as v — > oo. Requiring the absence of 
this divergence reduces the number of free parameters to 3. 

The above solutions Q( N \v) for N — 4 were used in Refs. [16], [33], [34] and [37] as pre- 
dictions for the scaling function of the order-parameter pair correlation function of three- 
dimensional systems at Lifshitz points. Furthermore, in Ref. [37] extensive Monte Carlo 
data were presented for the scaling function of the three-dimensional ANNNI model, 
which appeared to be in perfect agreement with these predictions. Since this case N = A 
is of particular interest to us, we give here the explicit form of the predicted Q( 4 '(v) for 
further use. It reads 

(2.20) 

(2.21) 
(2.22) 

(2.23) 

Aside from an overall (nonuniversal) amplitude b and the nonuniversal scale a, this 
scaling function involves a single universal parameter p = a 1 ' b\ /b . To adjust it 
by means of their Monte Carlo results for the three-dimensional ANNNI model, Pleim- 
ling and Henkel [37] considered ratios of truncated moment integrals Mj(vo) = J°° dv 
yi il^'foa 1 ' 4 ), where the use of a lower integration limit v > was necessary because 
they were unable to compute numerically the function Q(v) for values vo < 0.22. 

In the next section, we re-examine Henkel's arguments leading to the scaling func- 
tion (2.20). We will show that there are important reasons to question the presence of a 
contribution proportional to il\ (v) in il^'{v). 

3. Re-examination of the scaling-function solutions of the postulated 
invar iance equations 

Let us return to the postulated invariance equations (2.6)-(2.8). Unfortunately, it is 
not stated explicitly in Refs. [16,33,34] in what region of (t, r)-space these are presumed 
to hold. Clearly, in the case of a bulk equilibrium systems with a LP, the obvious point 
of view would be to interpret them as being valid in full d-dimensional space R x R d_1 , 
so that the t- variable is not restricted to positive values. 4 Accepting this interpretation, 
we can solve these equations by Fourier transformation. 5 Equations (2.7) and (2.6) yield 



rt 4 \v) = btf ) ntf\v) + b? ) n?\v) 

with 

o{ 4) r(3/4)^ IW2 + C/4) / -v 

l ' r(C/4)-^Z!r(Z/2 + 3/4) {2a 1 

and 

o(4) v^/2 ~ r[(* + i + o/4] *(0 

1 l ] r[(c + i)/4]^r(z/4 + i)r[(z + 3)/2] 

where s(l) is defined by 

s(l) = 2~ 1/2 [cos(Ztt/4) + sin(/7r/4)] cos(?tt / 


— V 


_(4a)V4_ 
/4). 



4 Obviously, this would be different for time-dependent phenomena such as relaxational processes where 
one must carefully distinguish between future and past time directions. 

5 Needless to say that the position-space functions G(t, r) can be trusted to belong to the space of 
tempered distributions (the dual of the Schwartz space <S(K d ) of rapidly decreasing C°° functions), and 
hence to have well-defined Fourier transforms. 



and 



pd k 



kd k 



2-k 



JV-l-JV-2, 



N 



p-d p + NA 



G(k,p) = 0, 



G(k,p) = 0, 



respectively, where the Fourier transform G(k, p) is defined by 

dk r d d - 1 

2^ J (2tt 
and A means the scaling exponent 



G(t,r) 



L P 

,d-l 



J(kt+p-r] 



G(k,p), 



A = l 

TV 



A 



1 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



The unique solution to the first-order partial differential equations (3.1) and (3.2) can 
be easily found by the method of characteristics. It reads 



G(k,p) = [p 2 +4aN- 2 (ik) 



N' 



-A 



(3.5) 



Note that in the special case of N = 4 and A = 1 , the result reduces to the usual form 
(p 2 + crfc 4 ) -1 of the free momentum-space propagator at a LP (with a = a/4) [see e.g. 
Refs. [38,39,27] and Eq. (5.1) below]. 

Let us assume that 4a(ik) N /N 2 > 0. For even N > 4, this is the case when a (— 1) N / 2 > 
0. The Fourier backtransform of the function (3.5) with respect to the variable p may be 
gleaned from Rcf. [44, p. 288]. Using this, one finds that the scaling function is given by 

l/N 



rt N \v)=C N ((,a) dn^-V/ 2 



' cos 



N 2 

4k*i 



A" 



with 



C N ((,a) 



7r -(d+l)/2 2 2-d+(C-l)/iV 



N 2 
4|a| 



(C-i)/JV 



l/N 



(* N/2 ) 



(3.6) 



(3.7) 



T[(d-l)/2-(C-l)/N] 
where K v {z) is the Macdonald function (modified Bessel function of the second kind). The 
function fiW(v) is a (particular) solution of the ordinary differential equation (2.11). By 
construction, it is the unique scaling function (up to scales) consistent with the validity 
of Eqs. (2.6), (2.7), and (2.10) in full (t,r)-space R d . 

Let us consider the case of our primary concern, N — 4 with a > 0, in more detail. 
The differential equation (2.11) then is of third order. Using Mathematica [45], one 
easily arrives at the three linearly independent solutions 



h(y) 
h(v) 



1-F2 
1-F2 
1-F2 



C 1 3 v 4 
4 

c 

4 

c 

4 



2' 4' 16a 

13 5 v 4 
+ 4' 4'4' 167 



2' 4' 2' 16a 



(3.8) 



Moreover, the power series of f2g (v) and Q[ (v) given by Eqs. (2.21) and (2.22) with 
N = 4 can be summed explicitly and the integral (3.6) for il ( - 4 - ) (w) be computed to obtain 
the results 



( i 



nr M =/. M - 2 ^>^> ). 



and 



^(v) = C 4 (C,a)2«- 3 )/ 4 r(5/4)r(C/4)^ 4) ( V ). (3.10) 

Note that the integrals (3.6) are even in v. Hence fl^ cannot have a contribution 
proportional to the odd solution ji. Since il\ has a term ex fa, it also cannot contribute 
to Q( 4 K That is, if the invariance equations (2.6)-(2.8) are taken to hold in full (t,r)- 
space M d , the coefficient b\ in Eq. (2.20) must vanish, as evidenced by our explicit 
result (3.10). 

There is a further reason by which a contribution oc tt\ (v) to Q,( 4 >(v) is ruled out: 

The functions fl (v) and fl[ (v) both satisfy the boundary condition (2.13) and hence 
vary as t>~^ as v — > +oo. Being even in v, the former also vanishes ~ 1 1^ | *= as v — > — oo. 
By contrast, £1^ (v) diverges exponentially in this limit. A convenient way to obtain its 
asymptotic behavior is to use its integral representation 

1 V ' (2 7 r) 1 /2r[(l + C)/4] LW ' J y ' 

with 

/•OO 

P(v) = / dKK (c " 1)/2 [e" m - cos(wk) - sin(vK)] i%_i) /4 (K 2 ). (3.12) 

Jo 
In the limit v — > — oo, we can replace the Bessel function K u—i) / 4 (k 2 ) by its limiting form 

■^/tt/(2k 2 ) e~ K , ignore the contributions from the parts of the integrand proportional 
to cos(|u|k) and sin(|v|K), and determine the dominant contribution of the remaining 
integral by expanding the argument of the exponential about the saddle point k = |v|/2. 
This gives 



P(v) « / dKK«-V/\ / JL- e «l-l-2« 2 



V/2 / dx (|^|/2)CC-3)/2 e v'/*-x' = 2 l-C/2^ | v |(C-3)/2 e « 2 /4. (3^3) 



It may be tempting to argue that the scaling form (2.10) should be modified by re- 
placing t by its absolute value in the scaling function i7^ 4 ', writing 

G(t,r)=r- 2A rt 4 \\t\r- 1 / 2 ). (3.14) 

In this way, the divergence of a contribution proportional to £l\_ (|i|r -1 ' 2 ) for t — *• — oo 

would be avoided. However, whenever the coefficient b\ of Cl\ does not vanish, such a 
correlation function fails to satisfy equation (2.7) in full (£, r)-space M. d . Indeed, applica- 
tion of the operator Y_\ to G(t, r) yields 

,(4) 

Y -i G ^)=-^£ m r5>{t) (3.15) 

rather than zero. Note that for reasons discussed in the paragraph following Eq. (2.8), 
this conclusion does not hinge on the definition of fractional derivatives d t ~ when N is 
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not a natural number. Using Hcnkcl's definition of fractional derivatives given in Ref. [34, 
Appendix A], one can in fact determine limjv_>.4 Y\-N/2G{t, r) in a straightforward fash- 
ion to recover the same result, Eq. (3.15). We are grateful to Malte Henkel (private 
correspondence) and one referee who both went explicitly through this analysis, confirm- 
ing our conclusion that the function (3.14) satisfies the inhomogeneous equation (3.15) 
rather than its homogeneous counterpart, Eq. (2.7) with N = 4. 

The proposed result (2.20) for the scaling function Q( 4 \ when re-interpreted according 
to Eq. (3.14), does therefore not satisfy the original equation given in Henkel's work 
[33,34]. Let us nevertheless accept the predictions for G(t,r) specified by Eqs. (2.20) 
(2.22) and (3.14) for the moment 6 and work out their consequences for the Fourier 
transform G(k,p). The results will be checked against explicit RG results in 4 — e dimen- 
sions and shown to be incompatible with the latter in subsequent sections. 

Rather than starting directly from these equations, it is more convenient to go back to 
the partial differential equations (2.6) and (3.15). Their Fourier transforms are Eqs. (3.2) 
(with N set to 4), and 



(pd k ~-k 3 d p )G(k,p) 



JA) 



aVr'A(l-2A)k-p 

P 



P -2A 



Here 



A = 



i d-l rr -(4A+l)/2 e -ie-r 



(3.16) 

(3.17) 



where e e R d l is a unit vector. Equation (3.2) yields the scaling form 

G(fc,p)=^ 2A . 9 (p 2 fc- 4 ). (3.18) 

Substituting this into Eq. (3.16) gives us a differential equation for the scaling function 
g(w), namely 

w 1/2 [-(a + Aw) g'(w) laAr 1 g(w)} = a b { f> A (1 - 2 A). 
This is solved in a straightforward fashion to obtain 



(3.19) 



g(w) 



a \ 
Aw J 

Aw J 



a 

4 



,, + -&< 4 > i 4(l-2A) 



du 



a 
Au 



A-l 



-A 



I^U-^-vW^-A^-A) 



(3.20) 

Here a is a constant which must depend linearly on the coefficients b and b\ of 
Eq. (2.20): a = ao&o + aib\ . The explicit expressions for the coefficients ao and a\ are 
not important for us. 

It is now easy to see that the contribution ex b\ is unacceptable. The function g(w) 
has the asymptotic behavior 

const w K - 2 A bf ] y/w. 



j(w) 



w-j-0 



(3.21) 



6 One motivation for this acceptance is the previously mentioned remarkably good agreement of the 
corresponding predictions for the order-parameter two-point function with the results of Monte Carlo 
simulations for the three-dimensional ANNNI model [37]. 
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Upon substituting this into Eq. (3.18), we arrive at the limiting form 

G(k,p) w const/c- 4A -2^6^ 4) p 1 - 2A fc- 2 . (3.22) 

But for nonvanishing momentum k the function G(k,p) must have a Taylor expansion 
in p of the form 

G(k>0,p) = G {k) + G 2 {k)p 2 + O( P i ), (3.23) 

where the behaviors Go(k) ~ fc~ 2A / e and G2(k) ~ fc-( 2 + 2A )/ e are dictated by scaling. 
To understand the behavior of G 2 {k), note that the scaling dimension of the descendant 
operator d r cf>(t,r) is A + 1. This translates into the stated behavior of G 2 (k) upon 
Fourier transformation. Furthermore, a fc-dependent second-moment correlation length 
£p(k) governing the decay of G(t,r) as r — > oo can be defined via its square 

vu p v ^n P =o 2 / d m te- lkt j d d - m rG{t,r) K ' 

where m = 1 in the uniaxial case we are concerned with. The length £,p(k) scales as k~ l ' e \ 
it has a finite value as long as k does not vanish. Thus, Gk (r) = j d m t e~ lkt G(t, r) must 
decay exponentially on the scale of £p(k) in the large- r limit, and all even moments 
f d d ~ m rr 2: > Gk(r) with j — 1,2,... should exist. This in turn means that G(k > 0,p) 
must be analytic in p 2 at p = 0. Analogous considerations show that G(k,p > 0) must 
be analytic in k at k = and expandable as 

G(k,p > 0) = g (p) +g 2 (p) k 2 + 0(fc 4 ), (3.25) 

where go(p) ~ p~ 2A and g 2 {p) ~ p~ 2A ~ 2e as p — >• 0. 

The contribution ex &]/ in Eq. (3.22) yields a term <~ p 1_2A fc -2 that is nonanalytic 
in p at p = 0. It is inconsistent with the expansion (3.23) and violates the mentioned 
analyticity in p. Consequently, this contribution must not appear, and the coefficient 6^ 
should be zero. 

4. Model, renormalization, and renormalization-group equations 

4.1. O(n) model for critical behavior at m-axial Lifshitz points 

According to Ref. [37], the Monte Carlo results obtained for the three-dimensional 
ANNNI model at its uniaxial LP could be fitted very well to the scaling functions pro- 
posed by Henkel for the case with TV = 4, i.e. 9 = 1/2. From the e-expansion results of 
Refs. [38,39] it is well known that the anisotropy exponents 8(d, n, m) of m-axial LP differ 
from 1/2 at order e 2 . One therefore expects that 9 ^ 1/2 also for the three-dimensional 
ANNNI model. However, the difference 9 — 1/2 appears to be small. Field-theory esti- 
mates based on the e expansion gave 9(3, 1, 1) ~ 0.47 [39,46], and Monte Carlo simulation 
data seem to be consistent with values 0.48 < 9 < 1/2 [42,37]. 

Our aim here is to determine appropriate two-point correlation functions by means 
of RG improved perturbation theory in d < d* (m) dimensions and compare the results 
with Henkel's predictions. To this end, we will consider the model used in the two-loop 
RG analysis of Refs. [38,39]. Its Hamiltonian is given by 
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H= d d x 



°\ (dl<t>) 2 + \(d,ct>y + 1 a^) 2 + i 2 + J i^ 



(4.1) 



where 0(x) = (^ a (x), a = 1, . . . , n) is an n-component order-parameter field. Pairs of a 
and /3 indices are to be summed from 1 torn and from to + 1 to d, respectively. Thus, 3 2 
denotes the Laplacian d 2 — d z ■ d z in 2-space, while (dp4>) 2 — {df>(t>a)dp4)a is the square 
of the gradient d r 4> in r-space. We shall investigate the two-point cumulants of the order 
parameter </> and the energy density <j) 2 (x)/2. They arc defined by 

5 ab G^°\x-x') = (Mx)M*')) cum (4-2) 

and 

, , , \ cum 

GtW{x-J) S (±<?{x)\4?&)) > (4-3) 

respectively. 

Despite the fact that 9 ^ 1/2, we shall be able to perform nontrivial checks of Henkel's 
predictions. The reason is the following. Building on the work in Refs. [38,39], we can use 
the e expansion about the upper critical dimension d* (to) to investigate the behaviors 
of G^ 0,2 ) and G^ 2,0 ) at the LP. The critical exponents these functions involve (such as 
6 and A) as well as their scaling functions have expansions in powers of e. Since the 
anisotropy exponent 9 starts to deviate from 1/2 not before second order in e, Henkel's 
predictions — if correct — ought to comply with the 0(e) results of the e expansion. In 
the case of the order-parameter cumulant G^ 2 '°\ the contribution of zeroth order in e of 
the scaling function is given by Landau theory; the leading corrections are of order e 2 
and encountered at two-loop order. By contrast, the one-loop term of the energy-density 
cumulant G^ 0,2 ' yields a contribution of zeroth order in e to its scaling function. 7 As wc 
shall see, the latter is inconsistent with Henkel's predictions. 

4.2. Renormalization 

To proceed, it will be necessary to recall some background on the RG analysis of the 
above model with Hamiltonian (4.1). The renormalization of the (dimensionally regular- 
ized) M-point cumulants (connected correlation functions) 

Giu l,a M (XI, • • • , X M ) = (4> ai (Xi) ■ • • Km (xm))™™ (4.4) 

of this theory in d = d* (to) — e < d* (to) dimensions has been explained in detail in 
Refs. [38,39,47]. Their UV divergences can be absorbed by means of the reparametriza- 
tions 

4> = Z l J 2 4> K , a = Z a o , u a~ m/i F m ^ = // Z u u , 
t~t lp = v 2 Z t (t + A tP 2 ), (p-plp) a- 1/2 = nZ p p, (4.5) 



7 An analogous situation is encountered in checks of the LSI predictions for type-II systems [36] . However, 
in this case the 0(e 2 ) term of the order-parameter response function could be determined explicitly and 
shown to be in conflict with the LSI prediction. Energy-density response and correlation functions were 
not computed in this work. By analogy with our case, the scaling functions of these dynamic cumulants 
will receive nontrivial corrections already at first order in e. The latter may be expected to show deviations 
from the LSI predictions as well. 
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where p is a momentum scale. Following these references, we choose the factor F m ^ as 

r(i + e /2)r 2 (i- £ /2)r(m/4) 

m ' e (4^)( 8 +™- 2e )/ 4 r(2-e)r(m/2)' l ' J 

The LP is located at (r,p) = (tlp,/3lp). In a theory regularized by means of a large- 
momentum cutoff A, the renormalization functions tlp and />lp would diverge ~ A 2 
and ~ A, respectively. However, in our perturbative approach based on dimensional 
rcgularization, they vanish. Results to order u 2 for the renormalization factors Zj,, Z a , 
Z p , Z T and Z u can be found in Eqs. (40)-(50) of Ref. [39]. The function A T is given to 
O(u) in Eq. (17) of Ref. [47]. Since we will not move away from the LP, we shall not need 
it. 

Using the reparametrizations (4.5), we can define renormalizcd M-point cumulants by 

C,..,^^-^.^.^^'''^....,,^.-;^.^). (4.7) 

(2 0) (2) 

When M = 2, this defines us the renormalizcd function G^' = G^ . However, the 
renormalization of the energy-density cumulant G^ ' 2 ' also involves an additive countert- 
erm. A possible way of fixing it is to subtract from & 12 ' its value at a normalization 
point (NP). We choose NP at the LP and a momentum (k,p) = (0,pp), where p is an 
arbitrary (d — m)-dimcnsional unit vector, defining 

G^ 2) (k, p; t, p, a, u, p) = Z 2 [& ^ (fc, p; f,p,a,u) - B(u,a; p)] (4.8) 

with 

B(u,v;p) = &°- 2 \ p ^G^ 2 \^ p ^ Lp . (4.9) 

Hence the renormalized function G^ ' satisfies the normalization condition 

G^ 2) \ Np = G^ 2 \o,pp;r = 0,p = 0,a,u,p)=0. (4.10) 

4.3. Renormalization- group equations 

The RG equations one obtains from Eqs. (4.5) and (4.7) upon varying p are known 
from Refs. [38] and [39]. Let us introduce the operator 

Vp=pd^+ J2 Wg (4.H) 

g=u,a,p,T 

along with the beta functions 

A; = MJo.9, g = u,<7,p,T, (4.12) 

where d^\o means a derivative at fixed bare interaction constant u and parameters a, p, 
and f . The j3 g can be expressed in terms of the exponent functions 

T) g (u) = pd^olnZg, g = (j),u,a,p,T, (4.13) 

and the function 

b T {u) = A T [pd^\ lnA r +T] T - 2r) p ] (4.14) 
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as 



u[e + r] u (u)}, 




g = u 


ar] a (u), 




g = v 


p[l + r]p(u)], 




g = p 


t[2 + 77 t (u)]- 


-b T (u)p 2 , 


g = r 



P 9 ={ I?., , ( 4 -!5) 



Two-loop results for the functions /3 U , 77^, r\ a , rj p , and r\ r may be gleaned from Eqs. (59), 
(58), and (40)-(50) of Ref. [39]. The function b T (u) (not needed in the following) is given 
to first order in u in Eq. (40) of Ref. [48]. 

The RG equations for the functions G R can be written as 



Vn + Y V*) G^ aN ; R ( Xl , ...,x N )=0. (4.16) 

Owing to the additive counterterm, the RG equation for Gp/ becomes inhomogeneous; 
it reads 

(D„ -2 Vt ) G^' 2) (k,p) = ^a-^Biu), (4.17) 

where B(u) is a UV finite function defined by 

p- f -a- m/i B(u) = -Z 2 T pd^B{u,a-p). (4.18) 

The RG Eqs. (4.16) and (4.17) can be solved in a standard fashion using characteristics 
[38,39,47,48]. For our purposes it will be sufficient to focus on the solutions to the RG 
equations of G R ' and G^ ' at the LP. The critical exponents they involve — namely, 
the correlation exponent 77L2, the correlation- length exponent z^L2, and the anisotropy 
exponent — may be expressed in terms of the values rj* = r] g (u*) of the exponent 
functions rj g (u) at the infrared-stable zero 

u* = ^^-+u* 2 e 2 + 0^) (4.19) 

3 n + 8 

of f3 u whose coefficient u?i, albeit not needed here, may be found in Eq. (60) of Ref. [39]. 
We have 

riL2=vl, ^L2-(2 + 77;) _1 , and0=(2 + 77;)/4. (4.20) 

Solving the RG equation for G R ' at the LP r = p = gives 

G^ 0) (fc,p) « M- 2 (£;)- 2 (p/A«) , ^- 2 * ro , d [(£:^) 1/ V 1/2 *(p//i)- tf ] (4-21) 

with 

*m,d(k) = G^ 0) (k,p; r = 0, p = 0, <r = 1, «*, a* = 1). (4.22) 

Here ££ = E* a (u) is a familiar nonuniversal amplitude. Just as its analog E*(u), which 

■y(0, 

r R 



we will encounter when solving the RG equation for Gk , it can be expressed as an 



integral along a RG trajectory: 

■Vg(u') 



E* (u) =E g (u*,u), Eg (u, u) = exp 



du' 



&(«') 



<t,t,4>. (4.23) 



The solution of the inhomogeneous RG equation (4.17) for the energy-density correla- 
tion function gives 

rto&fu „\~ (E*) 2 B(u*)v L2 f (p\~ aL/uL2 ^ [r „ t ,1/4 -1/2,/ / wen -,1 

(4.24) 
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with 

T m , d (k) = 1 + — ^ G { °' 2 \k,p;r = 0,p = 0,* = l,7i*, /i = 1) (4.25) 

o{u ) v h2 

where 

a L = 2- {d-m + m9)v L2 (4.26) 

is the specific heat exponent [39]. 

5. Two-loop calculation of scaling functions 

5.1. Calculations and results for general values of m 

We next turn to the calculation of the scaling functions (4.22) and (4.25) by means of 
RG improved perturbation theory. 

The Fourier transform of the free propagator Gf(x) reduces at the LP to the simple 
expression 

G f (k,p)= [p 2 +ak 4 y\ (5.1) 

Unfortunately, its Fourier backtransform yields for general values of m and d a rather 
complicated expression for the scaling function of Gf(x). One finds [38,39] 

G f (x) = r-^ a~ m ' A $ m , d (,T- 1 /4 zr -i/2) ( 5. 2) 

with 

Vm,d{V) 2 2 +m7r(6+rn _ 2e)/4 [r[( m + 2 )/4] 1 2 V 2' 2' 4 '64/ 

~ 4r(l + m/4) lF2 V^~ ; 2' 1 + T ; 64JJ- (5 ' 3) 

Owing to our use of dimensional regularization, the contribution from the one-loop 

graph v,./ vanishes at the LP. Hence the perturbation expansions to two-loop order of 
the cumulants & 2,0 \k,p) and G(°' 2 \k,p) at the LP become 

u 2 n + 2 
[W>-'"'{K,p)\ =p"+o-K^- 

and 



[GWfop)]- 1 =p 2 + ak i ^ H±± M k,p;a) + 0(u 3 ) (5.4) 



&°< 2 Hk,p) = \ J 2 (k,p;a) "" ( " 2 +2) [Mk,p;a)} 2 + 0(u 2 ), (5.5) 



respectively. Here Jj=2,3 are the integrals 

Jj (fe, p- Si) = f d d x [G f (x)} 1 e -i<P-r+k-*) 

_(a- m / 4 p-*I 2 (a 1 / A kp- 1 / 2 ), j = 2, 
" \a- m ! 2 p i -* e I s (a 1 ! 4 kp- 1 ' 2 ), j = 3. 



(5.6) 



Splitting off a factor F^ e from Ij(k; m, d) = Ij(k), let us write the Laurent expansions 
of the resulting ratios as 






RU[(k;m)e- 1 +R}? > (k;m) + RY , {k;m)e + 0(e 2 ), j = 2,3. (5.7) 
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In Appendix A we determine the low-order coefficients. As residues we recover the results 

of Refs. [38] and [39], namely 

R ( ?l(k;m) = l (5.8) 

and 

fl(3) (k;m) = J>)^ _ UW (5.9) 

- n ' ' 16to(to + 2) 2(8 - m) V ' 

Following Ref. [39] , we have introduced here the integrals 

r-OO 

j a (m) = B m / v m+3 ^ m4 ,{v)dv (5.10) 

Jo 



and 



with 



where 



U(m) = B m / t; m - 1 ^ ftI .(i;)dv (5.11) 

Jo 

R _ Srn-lS^-rn-l _ 2 10 + m 7r 6 + 3m / 4 r(m/2) 
m_ **„ " r(2-m/4)r»(m/4) ' ^ 



S m = 2^( m+1 )/ 2 /r[(m + l)/2] (5.13) 

is the volume of the unit TO-sphcre S m . Unfortunately, analytic results are known for the 
integrals j<j>(m) and Jo- (to) only for the special values to = 2 and 6, as well as for the 
limiting values ^(0+) and j CT (8-); see Eqs. (51)-(55), (84), and (86) of Ref. [39]. For 
general choices of to, these integrals can be determined by numerical integration. Results 
can be found in Table 1 of Ref. [39] . 

(2) 

In Appendix A we show that the regular parts i? ;>0 (k = 0; to) vanish exactly at k = 



i2, (2) (0;m) =*!,_!. (5.14) 



[Eq. (A. 10)]. Hence 

The finite parts R (k; to) and Rq (k; to) can be expressed as double integrals involving 
the generalized functions r^ and r^ , respectively. We have 

o6+m/2 m+4 roo 

4 2) (k;m) = r(2 + Tn/4)r(ro/4) y_ oo ^^Va.(r;k) (5.15) 



and 



R { 3) (k; to) = - i?^ (k; m) [2 7e - 4 + ip(2 - to/4) - ln(167r 3 )] 



B m / drr^ 3 A. (r;k), (5.16) 



where ip(x) = dlnT(a;)/da; is the digamma function [49]). Further, -je — — "0(1) is the 

) 

.,d* 



Euler-Mascheroni constant and f^ d * (r; k) means the function 



I'll , . n „ ( d — m r\ f°° , m i i , . „ /to rk i> \ 
^>5 k )= oFi ^^-,- T Jjf o ^^^^Mofi^y, —J (5.17) 

with <i set to d* = 4 + to/2. The reader may consult Refs. [44] and [18, Appendix] 
for general background on the distributions r^ and their Laurent expansions. The 
definitions of the generalized functions that are encountered in these expansions are 
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given in Eq. (A. 12) of Appendix A. Suffice it here to recall the simple example of r + . 
This distribution acts on test functions ip(r) as 

/OO /*1 />oo 

dr r^} <p(r) = drr~ l [(p(r) - tp(0)} + / drr~ X ip{r). (5.18) 

-oo JO Jl 

With the aid of the above results the renormalization functions B(u, <r; /i), B{u) 1 and 
the renormalized functions G^ ' and Gp ' can be computed in a straightforward man- 
ner. Our results are 

B(u,a;n)=fi a ' -— — — -5- + 0(u ), 

?(«) = -F m ,. 

>(3) 



B(u) = -F m , e - + 0(u), 
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R { '(k p ;m) 



-2 22^(fcp;m)ln(p//i) + 0(c) +0(u 3 ), 



(5.19) 
(5.20) 

(5.21) 



and 



4°' 2) (fe,P) = f-#^4{^ 2) fe;m)-ln^ + e4 2) (fe P 

£ Lb U ' I jLA 

-) 2i?[, 2) (fc p ;m)-ln. 



1/ r^ + 2^ 

— I € — 14 

2V 3 



«^[^ 2) (fcp;m)]^+O(r^.^.»0. 



(5.22) 

(5.23) 



where 

fc p = a^kp-V 2 . 

Upon setting u to its hxed-point value (4.19), one can convince oneself that the results 
are in conformity with the scaling forms (4.21) and (4.24). The scaling functions \& and 
T are found to have the e expansions 



l/* mid ._ e (k) = l + k 4 -2e 2 



n + 2 
(n + 8)'- 



R { o\k;m) + O(e 3 ) 



(5.24) 



and 



T m , d ,- £ (k) = l + ^(4 2) (k;m) + £i?f ) (k;m)-e^[4 2) (k;™)] 2 }+0(e 3 ) 
i/ L2 I n + 8 J 



? (2 



1 + J^ J (k; m) e + i?i j (k; to) e 



(2)/ 



«l/(^L2£) 



+ 0(e 3 ) 



In deriving the first and second form of Eq. (5.25), we made use of the fact that 



oil 4 — n 



e + 0(e 2 



(5.25) 



(5.26) 



^ L2 n + • 

The second form at this stage is just a convenient rewriting of the e expansion from 
which we shall benefit in Section 5.4 when extrapolating the scaling function T m =i,d*—e 
to d = 3 dimensions. 
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A cautionary remark is in order here. Expansions of scaling functions in powers of e 
such as those given in Eqs. (5.24) and the first line of Eq. (5.25), as well as those derived 
below, are not directly suitable for extrapolations to d = 3. They must be supplied with 
appropriate exponentiation hypotheses in order to capture the correct limiting behavior 
for k — >• oo. In the case of the functions 4 , m .d(k) and T TOi d(k), this asymptotic large- 
k behavior is dictated by the requirement that the p-depcndcncics of G R ' '(k,p) and 

G R (k,p) — G R '(0,/j,p) drop out. A convenient possibility to achieve the exponentia- 
tion of the asymptotic power laws in the limits k — > oo and k — > is to choose the scale 
parameter t = l\ such that the dimensionless inverse free susceptibility (p 2 + ak^/fj, 2 
becomes 1 at the chosen scale i\. Let us split off a factor p/fi, writing i\ — i^p/fj,. Then 
ly must be a solution to 

£^ 2 + £^ 4$ k 4 = 1. (5.27) 

Instead of Eq. (4.22), we then have 

tfrM(k) =tr- 2 G R 2 '°\e;%p^ 1 ;u*,l,l) (5.28) 

and a corresponding modification of Eq. (4.25), where 

(5.29) 

We shall return to the issue of the extrapolation of the scaling functions in Sec. 5.4 when 
extrapolating the e-expansion of T m=1 g/2- £ (k) to d = 3. 

5.2. The special cases d = m + 3, d = m + 3 = 5 — 2e, and (m, d) = (2,5 — e) 

For general values of m, the functions R\ appearing in the e expansions (5.24) and 
(5.25) of the scaling functions would have to be determined by numerical means. However, 
on the line d = m + 3, the functions & m .d( v ) reduce to simple Gaussians, 

$ m ,m+3(v) - (4^)-< m + 2 >/ 2 e- v2 l\ (5.30) 

This makes it possible to determine the functions i? ; (k; 2) in closed analytical form. 
Details of the calculations are described in Appendix A. To state the results and for 
subsequent use it is convenient to introduce the functions 

Q (s) =Im[(s-i)ln(s-i)] = -sti(s) - |ln(l + s 2 ) (5.31) 

with 

i?(s) = arctan(l/s) (5.32) 

and the dilogarithm Li2(z), a special one of the polylogarithm functions, defined by 
analytic continuation of their Taylor series (see Refs. [50] and [51, Sec. 2.6]). 

OO j 

Li fe (z) = ^^ for \z\ < 1 and k = 2, 3, . . . (5.33) 

3=1 J 

to the complex plane, with a branch cut along the positive real z-axis from z = 1 to 
z = oo. The analytical continuation of 1^2(2:) is provided by its integral representation 

[50-52] 
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Li 2 (*) = - I 

Jo 



dt 



ln(l - i) 
t 



1 , ln(l-zt) 
dt — 



We shall also need the Clausen integral (see Refs. [50], [52], and [53, Sec. 27. ■ 

sm(jz) 






r 



Im[Li 2 (e")] 



dt In 



2sm- 



(5.34) 



(5.35) 



The above special functions arc intensively used in papers dealing with calculations of 
Fcynman integrals appearing in usual <j> theory and QFT [54-58]. 
In terms of these quantities our results read 



?( 2 ) 



i?^(k;2) = Q (k72), 



(5.36) 



WM-S-i-fyrm-wm+Wn-** 



k 2 /4-l-k 4 

l[tf(k 2 /2)-tf(k 2 )]ln(i±^ 



>/W^ 



1 /4 + 4k 4 

2 ln UTk^ 



1 



■ln 2 (4 + k 4 ) + Rc 



Li 2 



2i + k 2 



jLi 2 (-k 4 /4) + |Li 2 (-k 4 ) - y Cl 2 [2^(k 2 /2)] 
+ ^ Cl 2 [2tf(k 2 )] + ^ Cl 2 [2tf(k 2 /2) - 2tf (k 2 )] , 



and 

i?f(k;2) 



18 V 6 



11 



IE Im 



3 



+ t 



9 



Im 



Y + t ) ln (ir + * 



(5.37) 



(5.38) 



Furthermore, it is possible to compute the required Feynman integrals on the line 
d = m + 3 in closed form (see Appendix A) for general values of m. When d = d* (m) — e, 
the constraint d = m + 3 implies that d = 5 — 2e and m = 2 — 2e. Using the results 
for the integrals Ij (k; m, m + 3) given in Eqs. (A. 40) and (A. 41), one can compute the e 
expansions of the scaling functions V&2-2e,5-2e and T 2 _ 2e .5_ 2e . Our results are 



l/*2- 2e ,5-2 e (k) = 1 + k 4 + 2e 2 -^U 4 [ 5 - In 11 



(n + : 
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Im 



k 2 V, a 2 . 

y + ln U +i 



64 
27 



3 



+ 0(e 3 ) 



(5.39) 



and 



T 2 - 2e , 5 - 2e (k) = l+6i-^{Qo(k 2 /2)+£Q 1 (k 2 /2)-£^i- [Q (k 2 /2)] 2 }+O(£ 3 ), (5.40) 
n + 8 I n + 8 J 



n + 8 
where Qi(s) is defined by 

7T 2 1 

Qi(«) = -g- + g Im t( s + *) ln ( s + *)] 

= ^ + 1 ln 2 (l + S 2 ) + £#(*) ln(l + ,s 2 ) - ^ 2 ( S ). 



(5.41) 
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In Fig. 1, the functions Q (k 2 /2) [Eq. (5.31)], i?^ 2) (k;2) [Eq. (5.37)], and Qi(k 2 /2) 

4r 




-1 - 



"755**^ 1.0 1.5 2.0 2.5 3.0 

"""^^ Qo(k 2 /2) 





Fig. 1. The functions Qo(k 2 /2) = R^\k;2) (blue, dashed), i^ 2) (k;2) (black, full line), and Qi(k 2 /2) 
(red, dotted). 

[Eq. (5.41)] are plotted as blue dashed, black full, and red dotted curves, respectively. 
Note that the 0(e) terms of the scaling functions Y2.5_ e (k) [Eq. (5.25)] and Y2-2e,5-2e(k) 
[Eq. (5.40)] agree and are given by eQ Q (k 2 /2) (4 - n)/(n + 8). The function Qi(k 2 /2) is 

(2) 

the analog of R{ (k; 2) for the case d— m = 3. Their difference is small. 
5.3. The energy- density cumulant in the uniaxial case to = 1 

Since the predictions of the LSI theory were used to fit the Monte Carlo data for 
the ANNNI model, the case of a uniaxial Lifshitz point to = 1 is of particular interest. 
Using our results for general values of m described in Section 5.1, we could determine the 
scaling functions for to = 1 by numerical integration. However, in order to compare with 
the predictions of the LSI theory it is preferable to have as much precise mathematical 
knowledge available as possible. It turns out that more detailed valuable analytical results 
can be obtained for the scaling function T of the energy-density cumulant. To do this we 
will start from the momentum-space representation of the Feynman integral /2(k; 1, d)> 
derive a contour integral representation of the latter, and relate it to solutions of a 
Fuchsian third order differential equation. 

The scaling properties of the integral J2(k,p;a) defined in Eq. (5.6) imply that the 
function /2(k; to = 1, d) can be written as 

/ 2 (k;l,d) = k- 2 V(k- 2 ) (5.42) 

with 

i- 00 dk f d d - 1 p 1 



J(P) J^ 2k J (2^-1 (pa + fc«)[(p - Pf + {k- I)*] ' ,r, - ,:$1 

The integral on the right-hand side of this equation converges for 5/2 < d < 9/2. It 
defines a function of P that is regular near the origin P = 0. Hence it has a Taylor 
expansion of the form 

oo 

J(P) = Y,CjP 2j - (5-44) 

3=0 
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On the other hand, the function P e J(P) is regular at P = oo. Therefore, J(P) can be 
expanded as 

oo 

J(P) = P^J2 B i p ~' (5-45) 

for large P. 

Much more information can be gained from the following contour-integral representa- 
tion of J(P) proved in Appendix C: 



c o4-2d 

Jrf-3 f 

(27r) d - 2 cos(?rd) 



J( p ) = TT^hi — 7T3T J( P ) (5-46) 



with 
where 

and 



J(w) - -i[Ji(w) + (1 - e- M ) J 3 (w)\ , (5.47) 

Ji(«;) = [ t+ " U) ^(t-l)^-^p(t,w) {d - 4)/2 (5-48) 

Jt_(iu) V^ 

J- 3 ( W ) = e 1 ^ 2 I' ^=(1 - i)( d - 4 )/ 2 p(t, ™) (d ~ 4)/2 . (5.49) 

Jo V* 

Here p(£, w) denotes the function 

/o(i,w) = (urf-l) 2 +4w (5.50) 

and 

t^{w)=w- 1 Tliw- 1/2 (5.51) 

are its zeros. 

The integration paths of the integrals (5.48) and (5.49) are illustrated in Fig. 2. For 
values of w with < w < 1, the path for Ji(w) is parallel to the imaginary axis. When 
1 < w < oo, the roots (5.51) have real parts Rei± = 1/w < 1. To define J\{w) in this 
case by proper analytic continuation, a path must be chosen that goes around the branch 
cut (0, 1) of its integrand. Deforming this path such that the subpaths along the Imt < 
and Im£ > rims of the branch cut pass through ±i0, respectively, one sees that the sum 
of integrals along these subpaths cancels the contribution to J{w) from Js{w). Hence 
we can rewrite S{w) as 

(I— iO nt + (w)\ J. 

/ +/ \- r [{t-l)p{t,w)] { - d - A ^ 2 when0<u;<oo. (5.52) 

Jt-(w) JiO J Vt 

The integrals (5.48) and (5.49) as well as (5.52) converge for d > 2. The UV singu- 
larities of J(P) therefore originate from the zeros of the factor cos(7r<i) = sin(7re) in the 
denominator of the right-hand side of Eq. (5.46). The integrals (5.48) and (5.49) are of 
the Euler type. They belong to the same class as Euler's hypergeometric integral 

r(c) f dw 

2 F 1 {a,b;c;z)= mr{c _ b) J o w i- b {1 _ w) i- c+b {1 _ wz)a ^ 53 > 

defining the hypergeometric function 2^1 [see e.g. Eq. (15.3.1) of Ref. [53]]. There are 
several ways to show that the so-defined function is a solution to the hypergeometric 
differential equation [59]. We pursue a similar strategy here. Following the procedure 
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lint 



t+(w<) 



t+(w>) 



t 



t-(w>) 



s- 



Ret 



t_( w< ) 

Fig. 2. Integration paths for the integrals J7i(ui) (black), J73(ui<) (blue), and S(w>) (blue) with 
< w < < 1 and 1 < u>> < oo, respectively. 

described in Masaaki's book [59, pp. 88-89], one can derive an ordinary differential equa- 
tion of third order and Fuchsian type [60,61] that is solved by the Euler integrals (5.48) 
and (5.49). Such differential equations are frequently exploited in studies of Feynman 
integrals and their singularities, see for example [62]. 

The functions J^(w), J7i(w;), and ^(w) all obey the same differential equation. Let us 
introduce the operator 

D = w — , (5.54) 

aw 

along with the parameter 

A=(d-4)/2 (5.55) 

and the coefficients 

A(l + 6A)(8A-9) A(3 + 4A)(1 + 6A) 



a = 



ai = 



6(1 + 4w) 6(1 + w) 

392A - 39 - 304A 2 A (52A - 5) (4A-3)(4A-1) 



48(1 + 4w) 



and 



2A 



6(1 + w) 

4A 



16 



a 2 



2(1+ Aw) 1+w 
Then this differential equation can be written as 



2A. 



(2 - 2A + D)(\ -2X + D)(D- 2\)J(w) 

= a Q J(w) +ai{D- 2\)J{w) + a 2 (1 - 2A + D)(D - 2\)J{w). 



(5.56) 

(5.57) 

(5.58) 
(5.59) 
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Inspection of the coefficients of its terms proportional to d k J/dw k with k = 3, 2, 1,0 
reveals that it is a Fuchsian differential equation with regular singular points at w = 0, 
w = —1/4, w = — 1, and u> = oo. It has three linearly independent solutions, two of 
which are regular at the origin. The pole of the coefficients ao and a\ that is closest to 
the origin is located at w = —1/4. Hence the Taylor expansions of solutions JJ(w) that 
are regular at the origin, 



J(w) 



3C 

^ AjW j , 



j=o 



(5.60) 



are guaranteed to converge inside the disc \w\ < 1/4. Substituting this expansion into 
Eq. (5.59) leads to the recursion relations 



A 



i+2 



(4A - 1 - 4j)(3 + 4j - 4A)(j - 2A) 
2(l+j)(2 + j)(5 + 2j + 2A) 
35 + 54 j + 20 j 2 - 26A - 12jA - 



A 



j 
8A 2 



2(2 + j)(5 + 2j + 2A) 



A 



j+i- 



(5.61) 



The low-order coefficients Aq and A\ can be computed in a straightforward fashion 
from Eqs. (5.47)-(5.49). The change of variables t — w~ x + 2iw~ 1 / 2 v — > v transforms 
J7i(w) into an integral from v = —1 to u = +1. Expanding the integrand in powers of 
yjw to order w and integrating term by term then gives 

- 2A + 8A 2 



JiH 



i 2 1+2A Ny, 



1 



with 

The expansion 



Js(w) 



2(3 + 2A) 

iV A = 7r 1/2 r(A + l)/r(A + 3/2) 
2A(5 + 4A) 



ok) 



A'kX 



N x 



1 



+ 0(w 2 ) 



3 + 2A 
can be obtained in a similar fashion. Writing the Taylor expansion of J^(w) as 

oo 

J(w) = N x J2 A J wi > 
we can combine these results with Eq. (5.47) to conclude that 



and 



Ai 



Ao = 2[4 A - 
2A - 8A 2 - 3 



i(ttA)] 



4A 



5 + 4A 



sin("7rA). 



(5.62) 

(5.63) 

(5.64) 

(5.65) 

(5.66) 

(5.67) 



Substituting expression (5. ( 
yields 



3 + 2A 3 + 2A 

for Ao along with Eqs. (5.63) and (5.55) into Eq. (5.46) 



Co = J(0) = 2 4 - 3 <V 3 - d >/ 2 



2 d + 16sin(7rd/2) 



(5.68) 



T[(d- l)/2] cos(?rd) 

for the coefficient Co appearing in Eq. (5.44). The result agrees with the value of 7(0, 1) 
given by Eq. (B.23) of Ref. [40] when m = 1. 

The coefficients Aj with j > 2 can be determined, on the one hand, from „4o an d 
A\ with the aid of the recursion relations (5.61). On the other hand, a general explicit 
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expression for them can be derived from Eq. (26) of Ref. [63], where the result of the 
inner p integration of the double integral (5.43) for J(P) is given in terms of the Appcll 
[64] function F\. The parameter e used in Ref. [63] is defined as e — 4 — D, where D is 
the dimension of the p-integral. Hence, we must set e = 4 — (d — 1) = 5 — d = 1 — 2A, 
p x = P and m.1,2 = (&T 1/2) 2 in the relevant Eq. (26) of this reference. Substituting the 
standard representation of the Appell function F\ as a double series, given in Eq. (20) of 
Ref. [63], into Eq. (5.43) one can integrate term by term over k to determine the series 
expansion (5.44) of the integral J(P). The series expansion (5.65) of J(P 2 ) then follows 
via Eq. (5.46). Our results for the coefficients are 

= 2sin(7rA) T(2 + 2A) (1/2 - 2A) 2j - ■ 

■*' A^ r(l/2 + 2A) (l-A),(3/2)/ ] 

x 3 F 2 (-X, 1/2,1/2 -2X + 2j;l-X+j, 3/2 + j;l), (5.69) 

where (c)j = T(c + j)/T(c) denotes the Pochhammer symbol. For j — and j = 1, the 
last expression reproduces Eqs. (5.66) and (5.67), respectively. Using Mathematica [45], 
one can check that these coefficients, Eq. (5.69), satisfy the recursion relations (5.61). 

From Eq. (5.69) one can obtain the e expansion of the coefficients Aj to 0(e) in a simple 
manner. Recalling that e = 1/2— 2A and taking into account that (e)2k = e(2/c— l)!+0(e 2 ) 
as e — > 0, one finds 

A, = 3\/2-\/2(7r/2 + 21n2)e + 0(e 2 ) (5.70) 



and 

Aj = e : 7T77r7T777r : (~1) J 3^(-l/4, 1/2, 2j; 3/4+j, 3/2+j; l) + 0(e 2 ), j > 1. (5.71) 



_ 3y2(2j 1). ( _^j _ P [_i /A -i /0 „-, q //!_,_„• q/o _l»-. i N i^ni',2\ 
' (3/4), (3/2), 



The e expansions of the first two coefficients _4j with j '• > 1 are given by 

^i = -^(4 + 7r-21n2)e + 0(e 2 ) (5.72) 

and 

-4 2 = ^[9 + 2(7r-21n2)]e + 0(e 2 ). (5.73) 

Higher-order terms of their e expansions can be obtained by means of the methods 
developed in Refs. [65] and [66]. 

Next, let us consider the asymptotic behavior of J^(w) for w — > oo. As is clear from 
Eqs. (5.42) and (5.46), the limiting form of J(w — > oo) is needed to determine the value 
of the function -?2(k; 1, d) at k = and check its consistency with the previously obtained 
result for /2(C); to, d) [39] recalled in Eq. (A. 10). The required information can be derived 
by analytic means directly from the integral representation given in Eqs. (5.46)-(5.49). 
To see this, note that Eq. (5.45) translates into a large-u; expansion of the form 

DC' 

j(w) = w a - i/4 J2 b j w ~ j/2 - ( 5 - 74 ) 

3=0 

To determine the asymptotic term oc Bo, we substitute the approximation t±(w) w 
±2zu> -1 / 2 into the integral (5.52), obtaining 
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(p2iw~ 1/2 r \ / a \ A 

e MA / +e- <wA / }dU- 1 / 2 (t 2 -- 

Jo J-2iw- l /2j V w 

= «; A - 1 / 4 2 2A+3 / 2 cos[7r(A-l/4)] f dxx' 1 ' 2 ^- x 2 ) x . (5.75) 

Jo 

The remaining integral can be performed. We thus arrive at the result 



*_^^_ 1/4)1 Ee™. 



(5.76) 



This can be combined with Eqs. (5.46) and (5.42) to compute /2(C); 1, d). The result agrees 
indeed with Eq. (A. 10). 

5.4. Extrapolation of the m — 1 energy- density scaling function to d = 3 

Before turning to a more detailed discussion of the above results, we will first use 
them to obtain an extrapolation of the energy-density scaling function Yi^k) to d = 3 
dimensions. Starting from the second form of Eq. (5.25), we combine Eqs. (5.42), (5.46), 
(5.74), and (5.76) to conclude that its term in square brackets can be written as 

VT 2f - J(k- 4 )/B = 1 + R {2 \k; 1) e + R (2) (k; 1) e 2 + 0(e 3 ). (5.77) 

Substitution of the left-hand side into Eq. (5.25) yields the asymptotic large-k behavior 
[ < -2q l /^ L2 _ -phe ex p 0ncn t m ay be recognized as the expansion of — a^/(p^20) to 0(e). 
Hence, the e-expansion result (5.25) is consistent with the expected fc-dependency of 
G R to this order. 

In order to cast our e-expansion result in a form that is well suited for extrapolating it 
to d = 3, we follow the strategy which led to the crossover scaling form (5.28) of G R . 

We exploit the behavior of G R ' (k,p; 0,0, it*, 1, 1) under scale transformation, choosing 
the scale parameter ly again as the solution to Eq. (5.27). The analog of Eq. (5.28) leads 
to the form 



2>iTi,3(I> 2 k) 



Jt k -4t,49-2\l QL >"'-' 



"k 



So k 2 < (1 ~ 26) 



(5.78) 



^0 * * k 

where D\ and D^ are two nonuniversal metric factors we introduced to adjust the am- 
plitude of T13 and the scale of the variable k. We have this freedom since we know from 

Eq. (4.24) that the scaling form of G R ' involves two such factors (related to E* and 
E*). We will fix them via the normalization conditions 

Ti, 3 (0) = 1 (5.79) 

and 

lim k aL/(,/L29) Ti. 3 (k) = 1. (5.80) 

k— >oo 

For Bo, we substitute its exact (d = 3)-value Ko|a=-i/2- Further, we set e = 3/2. Tak- 
ing into account that J"(0) = N\Ao, J(w) ~ Bow~~ f / 2 as w — > 00, and the limiting 
behaviors (5.29) of £k, one finds from the normalization conditions: 

26 



D 1 = l, 

D 2 = [N,A /B }- 2e/3 | A= _ 1/a - [2^r(3/4)/r(l/4)] ~ 26/3 . (5.81) 

For the critical exponents Eq. (5.78) involves, we use the d = 3 estimates a^jvi^.i — 
0.21 and 8 ~ 0.47 of Ref. [39]. The function J'(u') can be computed numerically from 
the integral representation 

rt+(w) 

J(w) = 2 1m dtt~ 1/2 (t- I)' 1 ' 2 p(t,w)- 1/2 for0<u><oo, (5.82) 

JiO 

to which Eqs. (5.47)— (5.51) simplify when d = 3. Likewise, we use numerical means to 
solve Eq. (5.27) for iy. The resulting scaling function Ti^k) one obtains in this fashion 
from Eq. (5.78) is depicted in Fig. 3 and compared with the scaling function TJ™?! d _ 3 (k) 
of a free massless field theory with action 

nt c m = lj-^j £p^ (p 2 + fc4 ) A V-*>- P ■ (5-83) 

The two-point cumulant one obtains from this action, 
l.Oh 




0.2 



2 4 6 8 10 12 14 

Fig. 3. Comparison of the extrapolated scaling function Ti d=3(k) (red full line) defined in Eq. (5.78) 
with the LSI prediction Tf|^(k) (blue dotted line) given in Eq. (5.84). For further explanations, see main 
text. 

(£k>, P >£k, P ) cam =P~ 2A T£S(fcp- fl *) (2n) d 5(k + k>) S(p + p<) , 

Ok) = (l + k*)- A , (5.84) 

agrees with what the LSI prediction for N = 4 (On — 1/2) and n = m = 1 given in 
Eqs. (3.14) and (3.20) yields upon normalization according to Eqs. (5.79) and (5.80) 
if the unacceptable contribution ex b x is dropped. The exponent 2A corresponds to 
«l/^L2- We therefore used the above-mentioned (d = 3, m = l)-estimate «l/^L2 — 0.21. 
As one sees from Fig. 3, the difference between our extrapolation Ti^ and Tf^ 3 is fairly 
small — the two functions differ by at most 2%. 

Finally, let us turn to a comparison of our e-expansion results for the scaling function 



"*~i 9/2-e with Henkel's prediction. Consider first the situation when b[ = 0. In this case 
the LSI prediction is given by T5 r g°, 2 _ £ . To obtain the e expansion of this function the re- 
sult (5.26) with n = 1 must be substituted for 2A. The O(e) contribution of this function 
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originates exclusively from the 0(e) term of this exponent. Hence the e expansion of this 
scaling function clearly differs at order e from our result given by Eqs. (5.25), (5.77), and 
(5.82). Next, consider the LSI prediction with b\'^0. Since the contribution of order 
e° is given by the Gaussian result T^ r g e , 2 , the coefficient b\ would have to be of order 
e. The additional 0(e) term of the LSI scaling function resulting from the contribution 
ex b x is incompatible with our e-expansion result. In fact, this incompatibility is not 
only quantitative but qualitative: The LSI scaling function g(w) given in Eq. (3.20) has 
a single nontrivial singular point located at w = 1/4. By contrast, the function J^(w), 
which our 0(e) contribution involves, was found to have additional singularities, namely 
branching points. 

Let us also show that our result for the energy-density correlation function G R ' (k, p) 
given in Eq. (4.24), unlike Henkel's b]_ ^ LSI result, has expansions in p 2 when k > 
and in k 2 when p > of the forms (3.23) and (3.25), respectively. To this end, we return 
to Eq. (4.24) and set /j, = a = E* = E* = 1 for notational convenience. Our result (4.24) 
then becomes 



G 



(0,2) 



, R '(k,p)^C lp - a ^^T m . d (k) + C 2 (5.85) 

with k = k/p , where Ci and C 2 are constants. Further, the function T m) <j(k) can be 
written as [cf. Eq. (5.78)] 






J(k 



-4j>4,0- 



Bo k 2e r r 



e(l-20) 



"L/(^L2e) 



, k = k/D 2 



(5.86) 



In the limit p — > at fixed k > 0, the momentum k — > 00 and iy varies as k ' according 
to Eq. (5.29). Therefore Eq. (5.86) reduces to 



T 



,a(k).« w 



J{k- 2 / 6 ] 

Bo k e / 



«L/(^L2e) 



(5.87) 



Substituting this into our result for the energy-density cumulant G R ' (k,p) yields 



G {0 ' 2) (k r,) ~ Cl 



j{p 2 k- 2 ' e D 2 2 /e ) 
Bo k*/ e D- e/e 



a-L/{v-L2i) 



+ c 2 . 



(5.88) 



As we know from Eq. (5.60), CT{w) has a Taylor expansion at w = 0. Using this we 
see that the right-hand side of (5.88) is analytic in p 2 at p — > and fixed k, and hence 
complies with the expansion (3.23). 

In the limit k — > at fixed p > 0, we have k — > and ty s=s 1 from Eq. (5.29). This 
implies 

J ( k_ 4) W(^) 



G (0 ' 2) (k r>) - — 



c 2 . 



(5.89) 



BoP e k 2 * 

We can now insert the large-w expansion (5.74) for Sf{w) to conclude that our result for 
G K ' (k,p) has indeed an expansion of the form (3.25). 
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6. Summary and discussion 

In this paper we reconsidered Henkel's LSI theory for type-I systems and performed 
careful checks of its predictions. A major motive for our work was the apparently very 
good agreement of the LSI predictions with Monte Carlo results reported in Ref. [37]. 
Our paper has two qualitatively distinct parts. The first part dealt with the consequences 
of the conjectured invariance equations (2.3)-(2.8). Accepting these equations as given, 
we reanalyzed their solutions. As we have shown in Sec. 3 and Appendix B for the cases 
N = 4 and N > 4, these equations generally have less physically acceptable solutions 
than anticipated by Henkel [16,33,34]. Specifically, in the case N = 4 that concerns the 
comparison with Monte Carlo simulation for the three-dimensional ANNNI model [37], 
the contribution from the second linearly independent function, fl x , must be discarded 
for the following reasons. If the scaling variable is taken to involve the coordinate differ- 
ence t\ —t<i rather than its absolute value, Q\ (v) diverges exponentially as v — > — oo and 
hence is unacceptable. We therefore made the replacement v — > \v\ considering the func- 
tion Q\ (|v|). Rather than being a solution to Henkel's original homogeneous equation, 
our Eq. (2.7), this function turned out to be a solution to the inhomogeneous Eq. (3.15), 
involving an inhomogeneity proportional to the derivative of the 8 distribution. More 
importantly, we found that the contribution ex b x entails a violation of general ana- 
lyticity requirements (as discussed at the end of Sec. 3). Hence it is unacceptable. Its 
omission [by setting the coefficient o{ = in Eq. (3.20)], on the other hand, implies that 
the scaling function of the momentum-space energy-density cumulant of the LSI theory 
reduces to that of a free theory with the action (5.83), namely the function T^°^ given 
in Eq. (5.84). 

In the second part of the paper we used RG-improved perturbation theory in 4+to/2 — e 
dimensions to determine scaling functions of the order-parameter and energy-density 
cumulants in momentum space to two-loop order. The results are given in Sec. 5. For the 
special choice d = m + 3, closed analytical expressions could be obtained for the scaling 
functions' series expansions to second order in e. For the case of primary interest, the 
uniaxial case m = 1, we managed to derive a countour-integral representation for the 
two-loop term of the momentum-space energy-density cumulant. 

We found that the predictions of the LSI theory generally are inconsistent with RG- 
improved perturbation theory in d = d* (m) — e dimension. Only at the level of Landau 
theory for the order-parameter cumulant and the one-loop approximation for the energy- 
density cumulant, where the LSI theory yields scaling functions of massless effective 
free-field theories, did we find it to be in conformity with our systematic expansions for 
proper choices of the exponents A and 6. However, as soon as we went beyond these orders 
to include nontrivial corrections to the scaling functions, the results did not comply with 
the LSI theory. In the uniaxial scalar (to = n = 1) case of the momentum-space energy- 
density cumulant we investigated in great detail, our e-expansion results for the scaling 
function turned out to be inconsistent with the LSI predictions irrespective of whether 
a contribution from the function fl[ is taken into account (b x ^ 0) or not {b\ =0). 

There are other observations concerning our two-loop results, which complement the 
evidence against the viability of the LSI theory provided by our e-expansion results. In 
Appendix D, we studied the behavior of the function Sf{w) that the two-loop contribution 
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to the energy-density cumulant involves in the complex w-plane. Unlike the LSI scaling 
function (3.20), which has a single nontrivial singular point (located at w = —a/4), the 
function JJ(w) was found to have additional singularities, namely branching points. The 
interested reader may find a detailed explanation of the branching behavior of this and 
the related functions S\{w) and 1/3 (w) in that appendix. Their behavior in the complex 
w-plane differs qualitatively from that of the LSI function. We admit that these functions 
merely appear in RG-improved perturbation theory. However, given their qualitatively 
different behavior in the complex plane, it seems highly unlikely to us that proper re- 
summations of the perturbation series might yield results in conformity with the LSI 
predictions, even if we did not know about the incompatibility with the e-expansion 
results. 

Finally, let us comment on the apparently excellent agreement of the LSI predictions 
with Monte Carlo simulation for the ANNNI model found in Ref. [37]. There are two 
problems with the LSI predictions used in the comparison: (i) they were based on the 
value 9 = 1/2, which may be a good approximation but differs from the RG estimate 
9 ~ 0.47 of Ref. [39] [and pretends that all corrections of order 0(e 2 ) and higher to 
this classical value sum to zero when e = 3/2]; (ii) a contribution proportional to the 
second linearly independent function, which we found to be problematic, was taken into 
account. Using a value of 9 different from, but close to, 1/2 in the scaling plot of the Monte 
Carlo results and dropping the non-free-field contribution to the LSI prediction will make 
the agreement presumably less striking, though it is not unlikely to remain reasonably 
good. If so, the situation would be reminiscent of the relatively good quality of the 
Ornstcin-Zcrnikc approximation for the order-parameter two-point cumulant at a bulk 
critical point in d — 3 dimensions whose reasons are twofold: the values of the correlation 
exponents rj{d = 3, n) are close to the classical one tjmf = and corrections to the zero- 
loop result for the scaling function are small. One important ingredient for the eventual 
good agreement of the LSI prediction with the Monte Carlo data is the small deviation 
of 9 from its classical value 1/2. It is conceivable that corrections to order-parameter 
scaling functions of the free LP theory are also small. In fact, our extrapolation for the 
scaling function Ti^ of the energy-density cumulant presented in Fig. 3 exhibits small 
deviations from the free-field theory LSI prediction T^g , which in turn agrees with the 
one- loop approximation for the scaling function T1.3. 

In summary, we can draw two important conclusions: First, LSI theory is definitely not 
valid in a mathematical precise sense in the checked nontrivial case of critical behavior 
at Lifshitz points. To our knowledge, the only cases in which its predictions are safely 
known to be exact are those trivial ones in which it reproduces the results of free massless 
field theories. Hence, its predictive power and viability appears to be rather limited. 
Second, the seemingly very good agreement with Monte Carlo data reported in Ref. [37] 
is probably due to the fact that corrections to the Ornstcin-Zcrnikc theory happen to be 
small in this case of a uniaxial LP in d = 3 dimensions. This need not be so in other 
cases. The good agreement may therefore be deceptive. 
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Appendix A. Calculation of the integrals .^(k) and i3(k) 

In this appendix we derive various results for the integrals Ij(k) = Ij(k;m,d) defined 
by Eq. (5.6), including their Laurent expansions. Combining Eqs. (5.6) and (5.2) yields 

Jj-(k) = Jj(k,p;l) = I d m Z f dd-m rr -j(2-e) $j ( zr -l/2) ^(p-r+k-z) ( A ^ 

We first perform the angular integrations in the subspaces M™ and R d_m . Let 

X(m;k) = — — [ dAe ik * (A.2) 

Sm-l Js™- 1 

be the average of the function exp(zk • z) with z = z/z over S m ~ l . The function X(m; k) 
may be found in the form of a Taylor series in Eq. (A. 17) of Rcf. [67]. This series can be 
summed to obtain the closed-form expression 

X(m; k) = F 1 (m/2, -k 2 /4) = (k/2) 1 -" 1 / 2 J m/2 ^(k) r(m/2), m > 0. (A.3) 

Using this result and making a change of variable z — > v = zr 1 ' 2 in one of the radial 
integrations, the integrals Ij can be written as 

Jj(k;m ' d) = ^ dr r3-«+0-D« v (J) {r ; k) (A.4) 

with 

A(r; k) = oFl (^, 4) f~ dv v^ *£,» oK (f , - T *f) , (A.5) 

where e = 4 + m/2 — d, as usual. 

Let us first determine the values of Ij at k = 0. When k = 0, the double integrals over 
r and v factorize. The r-integrals required for l2(0;m, d) and I^,{Q;m,d) are analytically 
computable; the results are 

« i r , x, •> , n 2 e - 1 r(e/2)r(2-m/4-e/2) ,. N 

dr^oF 1 [(d-m)/2,-r»/4] = V(2-m/4- e ) (A ' 6) 

and 

, 9, <t r , , w 9,1 2 2e ~ 3 r(e- 1) T(2 - m/4 - e/2) ,. N 

drr 2e - 3 F! rf-m /2,-r 2 /4 i- j l 7 7 ; , A.7 

/o I (3 — m/4 — 3e/2) 

respectively. The associated w-integrals are the special cases Joj(m, d) of the integrals 

/>oo 

J ld (m,d)= dvv m - 1+l & md (v), j = 2,3; 1 = 0,. ..,4, (A.8) 

Jo 

previously used in Ref. [38]. The first one, Jo, 2, is known from Eq. (82) of this reference 
for general values of (m, d): 

J (mrt- 2- 2 - e r 2 (l- e /2)r(2-m/4- £ )r(m/4) 

■Mm, d) - (2 ^r(2- £ ) ' (A ' 9) 
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The result can be substituted along with Eqs. (A. 6) and (5.13) into expression (A. 4) for 
l2(k;m,d) to obtain 

J 2 (0;m,d) = ^M, (A.10) 

in conformity with Eq. (38) of Rcf. [39]. 
In a similar manner one finds 

7 3 (0; m, d) = Jo, 3 (m, d) r(3d/2 _ m _ 3)r(m/2) • (A-H) 

Unlike Jo, 2(^1, d), the integral Jq^{tti, d) is not known in closed analytical form for general 
values of (m,d). However, it can be computed for given values of (m,d) by numerical 
means [38,39]. 

We now turn to the calculation of the Laurent expansion of the integrals Ij in e 
for general values of m. The right-hand side of Eq. (A. 4) may be read as the image 
(r+ 3 , ^m d) °f t nc r-dependent test function ^(r; k) under the map provided 

by the generalized function rf ' e 3 [44]. These distributions have the Laurent ex- 
pansions 

rt W ~ 23 = \ -! ( ^-4)! + r '~" + "+ 2J l ^ + + ° (£2) ' (A ' 12) 

where 6"' (r) is the l-th derivative of the (5-distribution and the other distributions are 
defined by (see e.g. Refs. [44] and [18, Appendix]) 

(r+ s In' r + , <p) = j™ dr r~ s (In r) 1 L(0) - ]T J ^f (0) 

(A.13) 



9=0 

r^ 1 d s ~V 



where 0(x) is the Heaviside step function. 

The action of <5^ 2j ~ A > on ^'.(r; k) can be computed in a straightforward manner. One 
obtains 

(<*> vSd(; k )) = ¥>£>! k) = Jo, 2 (m, d) (A.14) 

and 

/r// (3) , . u (3)//, n n Jj, 3 (m,d) 4 J . 3 (m,d) ,. K , 

Combining Eqs. (A. 4), (A. 6), (A. 12), and (A.14) then gives 

M^ = ^(rM + g m -ig,- m -i (r;lj ^ ( . k)) + 0(£) (A . 16) 

with 

^ / 2 ~ e r(2-m/4-c) ,. , w , 

c -' (m '^ r(i + e /2)r(L/L/2) (A ' 17) 

= 1 + ^ [7e - 2 In 2 - i>(2 - m/4)] e + 0(e 2 ). (A.18) 

In order to be consistent with Eq. (A. 10), the 0(e°) term of the first expression on the 
right-hand side of Eq. (A. 16) must cancel the k- independent contribution of the second 
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one. One easily checks that this is indeed the case and yields Eq. (5.36) for R '(k; 2), a 
function that vanishes at k = 0. 

Since the two-loop contribution to G^ 2,0 ) is quadratic in u, we split off a factor F^ e 
in J3. We thus arrive at the Laurent expansion 



,V(ro)k 4 J4>{m) 



e 



16to(to + 2) 2(8- to) _ 
^(2-m/4)-ln(16^ 3 ) ) , _ 3 (3) , 

1- o r + ( r + ' to k )) + °( £ )- ( A - 19 ) 



I 3 (k;m, d) 
F 2 



The terms in Eqs. (A. 16) and (A. 19) involving r^ 1 and r^ 3 , respectively, cannot 
straightforwardly be evaluated in closed form for general values of to. However, we know 
from Eq. (5.30) that the scaling functions $ m .d reduce to simple Gaussians on the line 
d = m + 3. This enables us to determine the functions </?„ d . in closed form for m = 2. 
One obtains 

and 

sinr 



J 3 ) f r - k ) - Smr e -k 2 r/3 (A21) 

^,5^k)- 6l44 ^ e . (A.21) 

To determine (r + " J , 932 5(5 k)), we compute (r + J ,</?2 5G k)), obtaining 



drr e -V^(r;k) = -^l- i r( e -l)Im[(k 2 /2 + z) 1 - e ] (A.22) 

and 

/■OO -1 

dr r 2 - Vg(r; k) = -——-T(2e 3) Im[(k 2 /3 + i) 1 " 6 ]. (A.23) 
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The results can be Laurent expanded in e. The expansion coefficient of the e° terms then 
give us the required quantities (r + ~ 3 , <p 2 5 (; k)) . This leads to the results for R (k; 2) 
and R (k; 2) given in Eqs. (5.36) and (5.38), respectively. 

To compute R{ (k; 2), we subtract from ^(k; to, d) its value at k = 0, defining 

J 2 (k; to, d) = I 2 (k; to, d) - J 2 (0; m, d). (A.24) 

The subtraction eliminates the pole term. The desired function follows from the Taylor 
expansion of -^(k; m,d) to 0(e); we have 

R[ 2 \k;m) = /2 (k;m) - [1 - 7S /2 + ln(2V^)]4 2) (k;TO), (A.25) 

where Rq (k;m) = ^(kjm, d*)/F m fi and i^ (k;m) = d f I 2 (k;m,d* — e)| _„ is the coeffi- 
cient of the 0(e) term of /2, while the coefficient in square brackets results from the 0(e) 
term of F m<e . 

We now specialize to the case to = 2 and start from 

7 2 (k;2,5-e) = (27r) (3 - e)/2 f dr r 3(£ " 1)/2 Ji_,(r) f d 2 v^ 2 25 _ e {v)(c lk - v ^ -l). (A.26) 
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For the Bessel function we substitute its e expansion. It is convenient to use the integral 
representation of the 0(e) term's coefficient —d v J v (r)\ v _ 1 , 2 /2 one obtains from (see e.g. 
[68, entry 2.3.5.3]) 



Mr) = 



(r/2)- 



I> + 1/2) ^ J_ x 



f dte trt (l-t 2 y- 1/2 , Rei/>-l/2, 



(A.27) 



by differentiating both sides and interchanging the integration and differentiation on the 
right-hand side. This gives 



J i - 1 (r) = 



1 



IE 



In- 



sinr 



— dz cos(rz) ln(l - z 2 ) + O 
2tt Jo 



(c 2 



(A.28) 



The last integral can be done analytically yielding the known result in terms of sine and 
cosine integrals [69, p. 74] quoted in Eq. (49) of [70]. 

The square of the scaling function $ 2 ,5-e(«) is treated in an analogous fashion. We 
replace $2,5-6 (v) by its expansion to 0(e) and use an appropriate integral representation 
for the term linear in e. A convenient starting point to find the latter is the integral 
represention of the scaling function 

o-m-l„(2e-6-m)/4 pi 

$ m , 4+ m/2- e H = tr,„ - , ^T77i / dtt 1 - £ (l-t 2 )( m - 6+2e >/ 4 e-*' 2 /4. (A.29) 



r[(m 



2e)/4] 



One way to prove it is to perform the integration with the aid of Mathematica [45] to 
obtain the explicit result (5.3) when Re(m + 2e) > 2 and Ree < 2 (which can then be an- 
alytically continued). Alternatively, one can Taylor expand the exponential in Eq. (A.29) 
and integrate term by term. The result is the Taylor series of & m ,d given by Eqs. (10) 
and (11) of Rcf. [39]. 

Setting m = 2 and computing the Laurent expansion of the integral in Eq. (A.29) to 
order e°, one can show that the scaling function satisfies 



$2,5-e(u) 



-2+e/2 



16r(l + e/2) 



-v>/4_V 



dte- tv2 / 4 \n(l-t 2 ) + 0(e 2 ) 



(A.30) 



To prove this, we rewrite the power (1 — £ 2 )~ 1+e / 2 f the integral's integrand as — dt(l — 
t 2 ) e < 2 /e and integrate by parts. The contribution from the boundary term can be rewrit- 
ten as the limit a —$■ 0+ of 

a" e 1 f 1 , 1 f 
= / dtt~ l ~ € . 

e e J a 

The integral on the right-hand side can be combined with one of the two integrals pro- 
duced via integration by parts to obtain 



dtt 



-l-e 



[(1 - t 2 ) e / 2 e' tv2 / 4 - 1] = 7£ + Ei(v 2 /4) - ln(« 2 /4) + O(e). (A.31) 



In the remaining integral (— w 2 /4e) L dt e~ tv / 4 (1 — t 2 ) c / 2 t~ e we expand the e-dependent 
powers to 0(e) and perform the two integrals that do not involve ln(l — t 2 ). Adding up 
all contributions and multiplying by the prefactor then gives the stated result (A.30). 
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The integral remaining in Eq. (A. 30) is analytically calculable and can be expressed in 
terms of the exponential integral functions Ei and Ei [53]. However, we found it more 
convenient to work with the integral representation (A. 30) for <I> 2] 5_ e , rather than with 
the analytic expressions in terms of special functions. Likewise, we prefer to use the 
integral representation of Jn_ e y 2 given in the first two lines of Eqs. (A. 28). 

Upon substituting them into Eq. (A. 26), the required Gaussian integrations over v can 
be done in a straightforward fashion, giving 



J 2 (k;2,5-e) =Qo(k2/2) r + e (7E + ln7r) 

-F2.0 L I 



dr s ^L( e -^/^l)lnr 



dzln(l — z ) dr 

2 ./o Jo 

f 1 , ln(l-z 2 ) 

2e / dz — K - -+ 

1 + z 2 



cos(rz) , _ rk 2 /2 



dr ■ 



r 

smr 



(-1 



r\c 



cxp 



1 + z 



+ 0(e 2 ) (A.32) 



where Qo(s) is the function defined in Eq. (5.31). Performing the remaining r integrations, 
one arrives at 



4 2) (k; 2) = Q 1 (k 2 /2) + -Wi(k 2 /2) + W 2 (k 2 /2) - g (k 2 /2) In 2, 
where Q\(s) is the function (5.41), while W\{s) and W 2 (s) denote the integrals 

Wi{s) = J dz ln(l - z 2 ) In (l + ^) 



and 



f 1 ln(l-z 2 ) (2s /l + : 

W 2 {s) = I dz ^— — rrr- s 4- arctan 



o 



(1 + z) 2 



'1 + z 



2s 



In 



1 



4s 2 



(1 + z) 2 



(A.33) 



(A.34) 



(A.35) 



To compute W±, we express ln(l — z 2 ) as a sum of logarithms ln(l±z). The integral can 
then be evaluated using Mathematica [45] . The result involves the dilogarithm function 
Li 2 [2i/(i + s)]. The real part of this expression can be written as Li 2 (2cos<p, ip) in the 
notation of Ref. [50], with ip — aig[2i / (i + s)} = arctans. According to its Eq. (5.17), it is 
given by Li 2 (2 cos (p,ip) — arctan 2 (l/s) = $ 2 (s) where 'd(s) is the function introduced in 
Eq. (5.32). To rewrite the imaginary part of the expression, we use the inversion formula 
for the dilogarithm (see e.g. Refs. [50, Eq. (1.10)] and [71, p. 652]), 



7T 2 1 2 

Li 2 (z) = — Li 2 (l/z) — -In z + 7rilnz. 



|arg(-z)| < 7T, 



(A.36) 



This gives 
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Im 



Lio 



2i 

i + s 



1, /1 + s 2 ,, 
= — In I — - — J [arctan(s) — ir\ + Im 



, 1 s 
Li2 '2 +2 2 



1 + s 2 1 

= — i?(s) In — h Cl 2 (2arctans) — -Cl 2 (4arctans), 



(A.37) 



where Eq. (5.5) of Ref. [50] was used to arrive at the second expression. Application of 
the duplication formula [50, Eq. (4.17)] to the two Clausen integrals 0^(0) along with 
the relation tt/2 — arctans = arctan(l/s) finally yields 



Wi(«) 



+ 4(1 - In 2) Q (s) - 2i? 2 (s) + 2s tf (s) ln[(l + s 2 )/4] - 2s Cl 2 [2i?(s)]. (A.38) 



The calculation of the integral W2 proceeds along similar lines but is more cumbersome 
and lengthier. Without entering into details, we just record our final result: 

W 2 {s) = -j ln\/2 In -^_£_ + \ n ^l + 4s 2 + (- - In 4) sd(s) (A.39) 



vT 



2i + 2s 



+ i Li 2 (-4s 2 ) 



+ Ms) - i?(2s)l s In " - ' " + Re Li 5 

L2s s - 

- jLi 2 (-s 2 ) + S - {Cl 2 [2i9(2s)] - Cl 2 [2tf(s)] + Ch[2d{s) - 2tf(2s)]} . 

From the above equations, the result for R 2 (k; 2) given in Eq. (5.37) follows in a straight- 
forward manner. 

It is also possible to compute the integrals 7 2 and ^3 on the line d = m + 3 for general 
values of m with < m < 2. One finds 



7 2 (k;m, m + 3) 



1 (k 2 - 2») 



to/2 



(k 2 + 2i) 



to/2 



and 



/3(k;?ri, m + 3) 



^to,i-to/2 1 - m/2 (1 - i) ro - (1 + t) T 

r(-m - 1) 



^25+2m 3l+3m/2 7r 2+m 



(k 2 - 3*) m+± - (k 2 + 3t) 



m+1 



(A.40) 



(A.41) 



Using these results and setting m = d — 3 = 2 — 2e, one can determine the scaling 
functions ^I' 2 _ 2e] 5_ 2e and T 2 _ 2e 5 _ 2e to 0(e 2 ) in a straightforward manner. The results 
are given in Eqs. (5.39)-(5.41). 

Appendix B. Integral representations for the scaling function Q{v) 

The differential equation (2.11) plays a central role in Henkel's theory. He expressed 
its general solution (2.14) in terms of the generalized hypergeometric functions 2-Fjv-i- 
However, it is quite difficult to analyze the behavior of the function QS N > (v) for v — > 00 in 
this representation. Here we derive an alternative integral representation for the general 
solution of Eq. (2.11), which is more convenient for this purpose. 

The main result of this appendix is the following theorem. 
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Theorem B.l Let N > 2 be an integer and £ > 1. Then the general solution of the 
differential equation 

/ d N ' 1 d \ 

Lf {v) ^^-^--Cv)f(v) = 0, (B.l) 

is given by a linear combination of the (N — 1) functions fi{v) defined by 

fi(v)=*ig(a l v)-g(v) (B.2) 

with 

g(v)= dkk vN ' 2 exp {N/2f/ N vk K v (k N/2 ), l = l,...,N-l, (B.3) 

Jo L J 

ai = cxp(2ml/N), 
and 

c-i 

" = — ■ 

Proof Upon integrating by parts and using the Bessel differential equation for the Mac- 
donald function, one can show that 

L [a t g{a x v)] = (2/N) 2 / N - 1 2" I> + 1). (B.4) 

This proves the theorem since the right-hand side in (B.4) is independent of I. 

Several remarks are in order here. Note, first, that the function g(v) defined by the 
integral in Eq. (B.3) is holomorphic in the complex w-plane. Second, if f(v) is a solution 
to Eq. (B.l), then the function tt(v) = f(va^ x ' N ), with a ^ an arbitrary complex 
number, solves Eq. (2.11). Third, since the functions fi(v a~ l ' N ) and v p F p provide 
two sets of TV — 1 linearly independent solutions of Eq. (2.11), they must be related 
by a w-independent matrix. This yields, on the one hand, integral representations for 
the generalized hypergeometric functions (2.15) and, on the other hand, expresses the 
integrals /; as linear combinations of the functions v v iFn-\- For given integer N — 
3,4,.. ., the coefficients of these linear combinations can be determined explicitly by 
integrating Eq. (B.3) using Mathematica [45]. Taking into account that 

K v (z) = \ £- e- z [1 + 0(1 /z)\ for |argz| <3tt/2, (B.5) 

z—>-co y 2z 

one can find from (B.3) the asymptotic behavior of the function g(v) as v — > oo via the 
saddle point method. Its limiting form depends on argw and N. It is governed either 
by a saddle point of the integrand of the integral in Eq. (B.3) or else by the integrand's 
behavior in the vicinity of the origin k = 0. Since the complete analysis of the asymptotic 
large- 1 v | behavior of the function g(v) for generic argw and N is rather involved, we will 
focus our attention on those special cases that are relevant for the v — > oo asymptotics 
of the function Q(v). 

Upon substituting the limiting large-z form (B.5) for the Bessel function of the in- 
tegrand of the integral (B.3), the integrand becomes F(k)exp[TZ(v,k)} with F(k) — 

^j2 k N(2„-l)/4 and 

K(v, k) = (N/2) 2/N vk - k N ' 2 (B.6) 

for large positive k. For given v > 0, TZ{v, k) takes its maximum on the integration path 
< k < +oo at the k- value 

k s (v) - (2/N) 2/N v 2 l^- 2 \ (B.7) 
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To obtain the asymptotic behavior of the function g(v) as v — ¥ oo, we can therefore 
expand lZ(v,k) about k s (v) to second order, replace F(k) by F(k s ), and extend the 
lower and upper integration limits of the resulting Gaussian integral to ±00. We thus 
arrive at the asymptotic behavior 



V ■ -^ DC 



L-C 



, N f 2 \ N l-Af + C (N-2 JV y .- 

•^) = M/zv^hv w ^ exp — - w «- [1 + 0( W -)]. (B.. 



N -2\N r l JV 



To compare this with Henkel's results, we set a\ = 1 in equations (4.23) and (4.24) 
of Ref. [34]. The comparison with Eq. (B.8) shows that the leading exponential large-u 
divergence given in Eq. (4.24) of Ref. [34] originates from g(v). Henkel's condition (4.25) 
[34], which we reproduce in Eq. (2.17) simply means that the function g{v) is required 
not to contribute to Cl(v) so that the leading contribution to Q,(v) as v — > 00 results from 
the term a\g{a\v) in Eq. (B.2). 

Consider next the large u-asymptotics of g{a\v) for integer values N > 4. To this end, 
we must study the behavior of the function (B.3) as v — > 00 at fixed argt> = 2tt/N. The 
asymptotic form can be derived along lines similar to those followed to obtain Eq. (B.8), 
provided the angle argw is small enough. For such argw > 0, the saddle point (B.7) 
is located in the upper complex fc-plane slightly above the real axis. The integration 
path must be deformed into the steepest-decent curve crossing the complex saddle point 
&s(ciM)- The resulting contribution may be gleaned from Eq. (B.8) by substituting 
v — ¥ <J\V. It reads 



N /2\(HCP, , i-iv+c 

x cxp f^_l v N/(N-2) e i2n/(N-2)\ £ + ^-1)) _ 

(B.9) 

and gives the asymptotic large-v behavior when N > 6. For TV = 6, the limiting form 
simplifies to 

g(a lV ) rj 7rV372 3- (1+c)/6 [wcxp(7r«/3)] (c - 5)/4 exp(2«w 3 / 2 /3). (B.10) 

Thus, g(uiv) decays indeed as ti(^ -5 )/ 4 in this case. By contrast, when N > 7, the 
right-hand side of Eq. (B.9) diverges exponentially as v — > +oo. In order to prevent such 
subleading divergences in the scaling functions 0(i>), further constraints must be imposed 
in addition to Eq. (2.17) on the coefficients b p ' in Eq. (2.14). Otherwise the algebraic 
decay assumed by Henkel does not apply. 

Turning to the case 4 < N < 6, we note that cos[2tt/(7V - 2)] < for N < 6. 
Therefore, the right-hand side of Eq. (B.9) decays exponentially and does not describe 
the asymptotic behavior as v — > oo. A straightforward analysis reveals that the leading 
contribution to the integral (B.3) results from vicinity of the end point k = of the 
integration path, giving 

g(a lV ) » -2^W|Y /2 ^+0(0 (B.ll) 
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for 4 < N < 6. In this case, all divergences of Q(v) are contained in g(v), which is in 
agreement with Hcnkel's numerical check [34]. 

Appendix C. Contour integral representation for one-loop function 

In this Appendix we show that the one-loop Fcynman integral J(P) introduced in 
Eq. (5.43) has the contour-integral representation given in Eqs. (5.46)-(5.49). Choosing 
Cartesian coordinates pi, . . . ,Pd-i such that P = P/P points along the pi-axis, we 
decompose p as, p — piP + p±. The resulting integrand of J(P) is a rational function of 
Pi having four simple poles located at ±iy/p±_ + fc 4 and P±i\Jp\ + (fc — l) 4 . To perform 
the integration over p\, we close the contour in the upper half plane. The result then 
becomes a sum of 2iri times the residues at the poles in the upper half plane Impi > 0. 
In the contribution from the residue with Re p\ = P, we make the change of variables 
k — 1 —¥ k. We thus see that J(P) can be written as 

J(P 2 ) = Re[Z(P 2 )] (C.l) 

with 

^P 2 ) = 77t4t7 i^ dk [d d -l pi _ ' 



X (fc - l) 4 - fc 4 + P 2 - 2iP (pi + /C 4 )V2 ' ( C ' 2 ) 

The next step is to split the integration over fc into two parts, one from — oo to 0, and 
a second one from to oo. This leads to 



with 



1(P 2 ) =2+(P 2 ) +2-(P 2 ) (C.3) 

dk f d d - 2 Pl _ 1 

; 4 )V2 

(C.4) 



2irJ (2ir) d - 2 (p^+fc 4 )i/ 
1 



(fc ± l) 4 - fc 4 + P 2 - 2iP{p\ + fc 4 )V2 ' 

where the positive values of both square roots are chosen. In the integral over p± we 
perform the angular integrations and make the change of variables pj_ — > 9 with pj_ = 
fc 2 sinh#. This gives 

C roc poo L,2d-6 

l ± (P 2 ) = 7^r d6 S mh d - 3 dk—— (C.5) 

(27r) d J J 9±{k,9) 

with 

g± (fc, 9) = (fc ± l) 4 - fc 4 + P 2 - 2iPk 2 cosh 9. (C.6) 

Let us hrst consider the integral Z_(P 2 ). To compute its inner integral in Eq. (C.5) 
by means of residue calculus, we combine the integrals along paths infinitcsimally above 
and below the real axis to conclude that 

foo pd-6 i f pd-6 

dk ^T7* = 1 — d / dk ^TTv (°- 7 ) 

io g-(k,V) l-e 4md J c g-{k,0) 

where C is the integration path shown in Fig. C.l. The integrand has three simple poles 

39 



Imk 

h(0) 



0.5 



-0.5 



G 



-0.5 -- 



-1.0 -L 



C 



Re A; 



0.5 



1.0 



• He) 



Fig. C.l. Integration path C of the integral on the right-hand side of Eq. (C.7) and poles kj{6) of the 
integrand. The integration contour C can be closed by a circle of radius R — > oo. 

away from the real axis, located at the zeros kj(9) of the cubic equation 



9-(k,0)\ k=k(0) = [1 - 4fc 3 - 2ik 2 Pcosh{9) + 6k 2 - 4fc + P 2 ] 
= 0, j = 1,2,3. 
We choose them in such a way that 
1 iP 



k=k D (6) 



(C.8) 



MO) 
Mo) 



fci(O) = 

i i 

2 + 2 
1 1 



7r < argfci(0) < 27r, 
(-1 - 2iP) 1/2 , n < arg k 2 {9) < 2tt, 



(C.9) 



y i (-\-2iP) l l 2 , 0<axgk 3 (9)<n, 

for < 9 < oo and real P > 0. 

Since the integration contour C can be closed by a circle of radius R — > oo, we can 
apply the residue theorem. Upon exploiting Eqs. (C.6) and (C.8), we find for the residues 



Res 

k=kj(9) 



1 



MM) 



S _(M) 



n -i 



dk 



= -k'A9) 



k=k j (6) 



dg-(k,0) 



09 



k=k 3 (9) 



It follows that 



2iP[k J (9)} 2 smht 



°° fc 2rf ' 6 _ ; 

„ g-(k,9) -(l-e^)P 



y» [M^ 8 



j'=i 



sinh# 



M*)> 



(CIO) 



(C.ll) 
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which inserted into Eq. (C.5), then yields 



3 



X - (P2) = 2P S (2^ 1^^X d9k ^ e) t^) sinh <" ( C - 12 ) 

The integral T + (P) can be dealt with in a similar manner. The poles of the /c-integral 
are now given by the zeros kj(6), j = 4, 5, 6, of the function g+(k, 9), namely 

k 4 (8) = 1^(6) e-™ , k 5 (6) = k 2 (0) e~ in , k 6 (0) = fc 3 (0) e iff . (C.13) 

Note that the chosen phases in Eqs. (C.9) and (C.13) guarantee that < argkj(9) < 2n 
for all j — 1, ..., 6 when < < oo and P > 0. The analog of Eq. (C.12) becomes 

C „ — 2i-nd 3 roc 

1 + (P 2 ) = - 2P{ t 3 )d -2 i^a^/ ^W[^)sinh<- 4 . (C.14) 

It is convenient to express the integrals over on the right-hand side of this equation 
in terms of the complex integration variables <;: 

<; = <Pj(6) = kj($) - 1/2, .7 = 1,2,3. (C.15) 

These maps tpj : (0,R) — > C parametrize paths Cj in the complex fc-plane. The total 
integral (C.3) therefore becomes an integral along the curve C3 — C\ — C2, 

I " , '° )f | ! ,XJ ,.«.*^" )l " (c ' 16) 

where 



P±(0 = \J- ^Iq P T (4-P 2 + 1 + 8^ 2 + 16c 4 ) (C.17) 

with argpj_(0) = 7r/2. As is illustrated in Fig. C.2, curves C\ and C3 start at ?i = 
fci(0) — 1/2 = — i/2 and ^3 = ^3(00) — 1/2, respectively, and terminate both at <j = —1/2. 
Curve C2 starts from c = &2(0) — 1/2 and runs towards 1 — ooi. Let us deform the path 
Ci into the union of paths C21 + C; m also displayed in Fig. C.2. Since &rgp±(<;) = on 
Ci m , the contribution from C; m to the integral J_ c d<; in Eq. (C.16) is purely imaginary. 
Thus it does not contribute to the real part of Z(P 2 ) and hence not to J(P 2 ) and J(P) 
[Eqs. (5.46) and (CI)]. The union of the remaining paths C3 — C\ — C21 can be deformed 
into a single path C32 drawn violet and dotted in Fig. C.2. Complex conjugation of the 
integral J„ gives an integral along the complex conjugate path C32, which starts at the 
complex conjugate ^3 of ^3 and terminates at c| . We thus arrive at the result 

Sd-3 f f _,_,__ /_m<2-4 1 „-iwd f . r / \ld-4 



w = iw^sll * b±wl + e k * M "H' (ai8) 

where again argpj_(0) = 7r/2. Finally, we transform from q to the integration variable 
t = -4<r 2 jw with w = P 2 . 

The integral representation (C.18) is equivalent to the one given by Eqs. (5.46)-(5.49). 
We shall prove this for values of P (~ 1) that are sufficiently small so that the absolute 
value r = \<&\ — \q\ satisfies r > |<ji|. The generalization to the half- axis < P < 00 
then follows by analytic continuation in P. 

Let us deform the integration path C32 of Fig. C.2 into the full violet curve depicted 
in Fig. C.3, and Cg 2 likewise into the dotted brown curve. Our choice of the phase of 
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Fig. C.2. Integration paths used in Eqs. (C.16) and (C.18). The location of the points ?i, q,. ..,<;% 
displayed in the figure corresponds to the choice P = 1. For further explanation, see main text. 



Im? 



•* If 

0.5OI . 



-1.0 



-0.5 



"3 



()., r ) 



Re? 



-0.5 i&?i 



-1.0 



Fig. C.3. Deformation of the integration paths C32 and C^ 2 into the full violet and dashed brown curves. 



p±(s) such that argpj_(ir — 0) = 7r/2 implies that the phase relations &rgp±(ir — 0) = 
argpj_(— ir + 0) = and aicgp±(ir — 0) = a,rgp±(—ir + 0) = hold on the paths C32 and 
C32, respectively. As a consequence, the sum of the contributions from the two integrals 
between ir =F and cj =p cancel. Likewise, the integrals along the paths between ?i =p 
and — ir =p add up to zero. The contributions from the remaining portions of the paths 
in Fig. C.2 add up to 



42 



j(p)= onf „ fe-7 J r*KW] d - 4 + r*KW] d - 4 



8P(27r) d - 2 cos(Trd) 

+ (l-e-- d ) r*[px(?)] d -n ) (C.19) 



where phases of p_i_(?) in the three integrals are fixed by the conditions 

a,rgp±(ir) = a,igp±(—ir) = 0, argpj_(0) = ir/2. (C.20) 

Further, the principal value of the power [pj_(<r)] d_4 is to be taken, i.e. arg[pj_(cr)] d ~ 4 = 

(d — 4) argpj_(cr). Since the integrands of these integrals are even, Eq. (C.19) can be 

rewritten as 

c ( rS2 rsi 

°d—3 I / , r / Md-4 , r-\ „-iird\ / » r /_\id-4 



4P(27r) d cos(Trd) [ 7 ? » 7 J 

Both integration paths lie in the lower half-plane Im? < 0, and the phase conditions 
(C.20) must be taken into account. Changing to the integration variable t — — 4^ 2 /w 
with w = P 2 finally yields the representation given in Eqs. (5.46)— (5.49). 

Appendix D. Branching of complex integrals Jk(w) 

The functions J\{w), ^(w), and J(w) we considered in Sec. 5.3 for < w < oo can 
be analytically continued to the complex w-plane. These analytic continuations become 
multivalued functions with four branch points at w\ — — 1, w 2 = — 1/4, W3 = 0, and 
u>4 = oo. The branching of these functions, which we are now going to study, is described 
by the monodromy group. 

To this end, consider a real value wq of w with —1/4 < wo < and arg wo = n. 
For such Wq, all branch points of the integrands of the integrals in terms of which the 
functions J\{w), J%, and J{w) were expressed in Eqs. (5.47)— (5.49) are real and given 

by 

ti(wo) = t-(wo), t 2 (w ) = t+(w ) , *3 = 0, <4 = 1, (D.l) 

where tzp(w) are the zeros (5.51) of the function p(t,w) introduced in Eq. (5.50). They 
satisfy t 1 (w ) < t 2 (w n ) < i 3 < t 4 . 

Let J7fe(w_), with k = 1, 2, 3, denote the integrals 

J k (w-) = w 2 _ x /* " fl rfi[t-ti( W _)] A [t-t 2 ( w _)] A (i-t 3 ) _1/2 (i-* 4 ) A . (D.2) 

Jt k 

For k — 1 and 3, these definitions comply with Eqs. (5.48) and (5.49). The integrand 

f-V 2 FJ fc=1 2 4 (t — tk) x is a multivalued function. Following Ref. [59], we fix its phase by 

setting 

(0 if j<k, 

arg(i -*,-)=< (D.3) 

[ -7T if fc + 1 < j. 

This guarantees that the integrand is an analytical function in the lower half-plane 
Imi < 0. 

The functions J7fc(u>_) can be analytically continued from the interval (—1/4,0) into 
the complex w-plane punctured at the branch points w% = —1,-1/4,0. At these branch 
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points, some of the endpoints tk{w) of the integration paths in Eq. (D.2) merge or become 
infinite. Namely, to,{w) —¥ £3 as w — > —1/4, ti t 2(w) — > oo as w — > 0, and t\(w) — > £3 as 
w — > — 1. To study the branching of the integrals J^k(w) at the points w\, W2, and 
W3, let us change w continuously by moving along loops 7$ that emanate from w- and 
terminate there, passing counter-clockwise around one of the branch points Wi of the 
functions Jk{w), as is illustrated in Fig. D.l. 

\raw 




\ — ►Rew 



0.50 -L 



Fig. D.l. The branching of the integrals j7fc(u>o), k = 1,2,3, is characterized completely by the mon- 
odromy transformations A (7;) pertaining to the loops 71, 72, and 73. 

As w is changed continuously, the functions ,7fc(u>) also change continuously. However, 
because of the nontrivial monodromy, they do not normally return to the starting values 
J(w-) if the loop is traversed a single time. Let A(ji) Jk{w_) denote the end value one 
reaches from J7fc(w_) by going once around the loop 7,. Although A(ji) Jfk{w~) generally 
differs from Jk (w— ) , it must be a linear combination of the three linearly independent 
solutions of the differential equation (5.59) it solves. One finds 
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with 

c A = e 2 " A and b x = (c x - l)(c 2 x + 1). (D.7) 

These equations describe the action of the monodromy group on the three basic inte- 
grals J7fe, k = 1, 2, 3. As an immediate consequence, we obtain for the monodromy group 
action on the integral J(w-) the result 

A(j 1 )J{w-) = J(w-)+ib x cx[Mw-)-J 3 (w-)}, 

A(j 2 )J{w-) = J(w-)-i[c x + c x - 1 ]j 2 {w-), 

A(k)J(w_) = J(w-). (D.8) 

Directly at the upper critical dimension d*(l) = 9/2, one has A = (9/2 — 4)/2 = 1/4, 
c x = i, and 6a = 0, as a consequence of which A( , -fk)J'(w-) = J(w-) for k = 1, 2, 3. 
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